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Abstract 
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This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS)  for  the^ period, October  1989  -  31  March  1990.  It  also 
presents  statistics  from  operation  of  the  Intelligent  Monitoring  System  (IMS). 
Under  an  agreement  between  NORSAR  and  the  University  of  Helsinki,1 £<}ata 
from  the  Finnish  Experimental  Seismic  Array  (FINESA)  are^being  recorded 
and  processed  at  the  NORSAR  Data  Center  as  well  as  in  Helsinki,  and  results 
from  the  NORSAR  processing  of  these  data  are  reported. 

The  NORSAR  ^Detection  .Processing  system  has  %een 'operated  through¬ 
out  the  period  with  an  average  uptime  of  91.8  as  compared  to  94.9  for  the 
previous  reporting  period.  A  total  of  1648  seismic  events  have  been^reported 
in  the  NORSAR  monthly  seismic  bulletin.  The  performance  of  the  continu¬ 
ous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC  has 
been  satisfactory.  Processing  of  requests  for  full  NORSAR/NORESS  data  on 
magnetic  tapes  has  progressed  according  to  established  schedules.  >£ - 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data 
Processing  Center  (NDPC)  of  NORESS,  ARCESS  and  FINESA  data  have 
been  conducted  throughout  the  period,  with  an  average  uptime  of  96.9%  for 
NORESS,  95.3%  for  ARCESS  and  98.5%  for  FINESA.  The  Intelligent  Mon¬ 
itoring  System  was  installed  at  NORSAR  in  December  1989  and  has  been 
operated  experimentally  since  1  January  1990.  Preliminary  results  of  the  IMS 
analysis  are  given. 

There  have  been  no  modifications  made  to  the  NORSAR  data  acquisition 
system.  The  process  of  evaluating  technical  options  for  upgrading  the  array  is 
continuing.  A  test  of  a  full  subarray  acquisition  system  will  be  performed  at 
the  end  of  the  next  reporting  period. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  SUN-3/280  data  acquisition  sys¬ 
tems.  The  routine  processing  of  FINESA  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  modification  of 
equipment  for  FINESA  and  preparatory  work  in  connection  with  NORESS 
HF  instrumentation.  Other  activities  involved  testing  of  the  NORSAR  com¬ 
munications  systems. 

During  14-17  February  1990  NORSAR  hosted  an  international  symposium 


ii 


entitled  “Regional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification”. 

The  symposium  demonstrated  the  considerable  progress  in  the  field  of  seis¬ 
mic  monitoring  during  recent  years.  It  particularly  highlighted  the  technolog¬ 
ical  advances  in  seismic  instrumentation,  data  communication  and  computer 
processing,  as  exemplified  by  the  development  of  advanced  regional  seismic 
arrays  with  very  sophisticated  automatic  and  interactive  signal  processing  fa¬ 
cilities.  The  presentations  at  the  scientific  symposium  show  that  these  techno¬ 
logical  advances  are  accompanied  by  considerable  scientific  progress,  although 
much  work  remains  in  order  to  fully  exploit  the  potential  offered  by  regional 
arrays  in  a  seismic  monitoring  context.  The  majority  of  the  papers  have  been 
submitted  for  publication  to  the  Bulletin  of  the  Seismological  Society  of  Amer¬ 
ica ,  and  will  be  published  in  a  Special  Issue  of  the  Bulletin,  scheduled  to  appear 
in  the  fall  of  1990. 

We  have  continued  our  work  aimed  at  evaluating  the  stability  of  RMS  Lg 
for  yield  estimation  purposes.  Data  for  Semipalatinsk  explosions  from  the 
British  station  GAM  in  Garm,  USSR,  have  been  analyzed  and  compared  to 
previously  reported  recordings  at  NORSAR,  IRIS  and  CDSN  stations.  These 
studies  confirm  that  for  explosions  at  Semipalatinsk  with  good  signal-to-noise 
ratio  (>1.5),  mfcfLg)  may  be  estimated  at  single  stations  with  an  accuracy 
(one  standard  deviation)  of  about  0.03  magnitude  unit.  It  is  noteworthy  that 
this  accuracy  is  consistently  obtained  for  a  variety  of  stations  at  very  different 
azimuths  and  distances,  even  though  the  basic  parameters  remain  exactly  as 
originally  proposed  by  Ringdal  for  NORSAR  recordings  (0.6-3.0  Hz  bandpass 
filter,  RMS  window  length  of  2  minutes,  centered  at  a  time  corresponding  to 
a  group  velocity  of  3.5  km/s).  In  particular,  comparison  of  GAM  and  ARU 
data  suggests  that  this  consistency  applies  over  a  range  of  two  full  magnitude 
units. 

A  study  is  presented  on  integrated  wavefield  decomposition  using  three- 
component  seismograms  and  array  data.  The  array  beam  parameters  are  used 
to  produce  estimates  of  the  P,  SV  and  SH  contributions  as  a  function  of  time. 
This  decomposition  of  the  seismic  wavefield  by  wavetype  as  a  function  of  time 
not  only  has  considerable  benefits  for  the  recognition  of  seismic  phases,  but 
also  provides  a  domain  in  which  the  relative  proportions  of  P,  SV  and  SH 
can  be  compared  directly,  because  free-surface  amplification  effects  have  been 
removed.  This  information  on  the  current  proportions  of  different  wavetypes 
summarizes  much  of  the  propagation  processes  between  source  and  receiver 
and  therefore  can  be  beneficial  in  attempts  to  discriminate  between  different 
source  types. 

A  study  on  interpreting  Lg  azimuth  anomalies  observed  at  NORESS  in 
terms  of  ray  calculations  is  reported.  Initial  attempts  to  relate  ray-calculated 
and  observed  azimuth  anomalies  have  yielded  somewhat  inconsistent  results, 
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but  upon  careful  examination  of  the  ray  behavior  it  becomes  evident  that 
the  nature  of  the  Moho  heterogeneity  in  the  vicinity  of  NORESS  is  probably 
largely  responsible.  The  complexity  of  the  Moho  saddle  east  of  NORESS,  its 
comparatively  complicated  geometry  over  a  region  of  restricted  spatial  extent, 
results  i n  ray  diagrams  which  are  very  sensitive  to  source  location  and  phase 
velocity.  The  study  suggests  that  the  Moho  saddle,  a  feature  which  appears  to 
be  well  defined  and  documented  in  a  number  of  independent  investigations,  will 
play  a  significant  role  in  distorting  the  wavefield  and  is  probably  responsible 
in  large  part  for  the  Lg  azimuth  anomalies  observed  at  NORESS  from  sources 
to  the  east. 

A  regionalized  time-of-day  analysis  of  six  months  of  IAS  data  from  the 
regional  arrays  NORESS  and  ARCESS  has  shown  that  about  97%  of  the  more 
well-located  events  probably  are  of  man-made  origin.  Based  on  the  derived 
information  about  where  and  when  the  man-made  events  occur,  contour  maps 
(in  magnitude)  of  associated  capabilities  for  reporting  natural  earthquakes 
are  provided,  for  different  times  or  hourly  intervals  of  the  day.  The  magnitude 
limits  vary  from  about  3.0  as  a  maximum  in  some  mining  areas  and  down  to 
less  than  1.0  for  western  and  northern  Norway  offshore  areas,  where  most  of 
the  man-made  disturbances  are  still  below  the  detection  level. 

The  concept  of  threshold  monitoring,  introduced  by  Ringdal  and  Kvaerna, 
is  a  method  of  monitoring  the  seismic  amplitude  levels  for  the  purpose  of  using 
this  information  to  assess  the  largest  size  of  events  that  might  go  undetected  by 
a  given  network.  In  an  effort  to  demonstrate  the  capabilities  of  this  threshold 
monitoring  concept,  a  preliminary  version  has  been  implemented  into  the  Intel¬ 
ligent  Monitoring  System  (IMS).  Using  recordings  by  the  NORESS,  ARCESS 
and  FINESA  arrays,  the  method  has  been  applied  for  monitoring  two  selected 
target  regions  for  a  24-hour  interval:  a  mining  site  near  the  Finland-USSR 
border  and  the  nuclear  test  site  at  Semipalatinsk.  It  is  demonstrated  that 
the  implementation  of  the  threshold  monitoring  method  in  the  IMS  system 
enables  real-time  operation.  The  displays  provided  by  the  threshold  monitor 
appear  to  be  very  valuable  in  pointing  out  time  intervals  of  particular  interest, 
thus  aiding  the  analyst  in  his  work.  The  interesting  intervals  can  be  examined 
by  different  processing  techniques  to  locate  and  identify  the  events.  The  ex¬ 
amples  in  this  study  have  demonstrated  that  the  method  can  be  applied  both 
at  regional  and  teleseismic  distances.  We  note,  however,  that  some  additional 
research  needs  to  be  done  to  assess  the  potential  bias  in  the  upper  magnitude 
limits  when  detected  phases  occur  from  events  in  the  target  region. 
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1  Summary 


This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic  Ar¬ 
ray  (ARCESS)  for  the  period  1  October  1989  -  31  March  1990.  It  also  presents 
statistics  results  from  operation  of  the  Intelligent  Monitoring  System  (IMS). 
Under  an  agreement  between  NORSAR  and  the  University  of  Helsinki,  data 
from  the  Finnish  Experimental  Seismic  Array  (FINESA)  are  being  recorded 
and  processed  at  the  NORSAR  Data  Center  as  well  as  in  Helsinki,  and  results 
from  the  NORSAR  processing  of  these  data  are  reported. 

The  NORSAR  Detection  Processing  system  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  91.8  as  compared  to  94.9  for  the 
previous  reporting  period.  A  total  of  1648  seismic  events  have  been  reported 
in  the  NORSAR  monthly  seismic  bulletin.  The  performance  of  the  continu¬ 
ous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC  has 
been  satisfactory.  Processing  of  requests  for  full  NORSAR/NORESS  data  on 
magnetic  tapes  has  progressed  according  to  established  schedules. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data 
Center  (NDPC)  of  NORESS,  ARCESS  and  FINESA  data  have  been  conducted 
throughout  the  period,  with  an  average  uptime  of  96.9%  for  NORESS,  95.3% 
for  ARCESS  and  98.5%  for  FINESA.  The  Intelligent  Monitoring  System  was 
installed  at  NORSAR  in  December  1989  and  has  been  operated  experimentally 
since  1  January  1990.  Preliminary  results  of  the  IMS  analysis  are  given. 

There  have  been  no  modifications  made  to  the  NORSAR  data  acquisition 
system.  The  process  of  evaluating  technical  options  for  upgrading  the  array  is 
continuing.  A  test  of  a  full  subarray  acquisition  system  will  be  performed  at 
the  end  of  the  next  reporting  period. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  SUN-3/280  data  acquisition  sys¬ 
tems.  The  routine  processing  of  FINESA  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  modification  of 
equipment  for  FINESA  and  preparatory  work  in  connection  with  NORESS 
HF  instrumentation.  Other  activities  involved  testing  of  the  NORSAR  com¬ 
munications  systems. 

The  research  activity  is  summarized  in  Section  7.  Section  7.1  presents  a 
report  from  the  symposium  on  “Regional  Seismic  Arrays  and  Nuclear  Test  Ban 
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Verification”  which  was  held  in  Oslo,  Norway,  14-17  February.  An  analysis 
of  data  from  the  British  station  GAM  near  Garm,  USSR,  for  Soviet  nuclear 
explosions  is  presented  in  Section  7.2.  A  study  of  wavefield  decomposition  for 
three-component  seismograms  is  given  in  Section  7.3.  Section  7.4  discusses 
a  ray-based  interpretation  of  Lg  azimuth  anomalies  at  NORESS.  A  study  of 
earthquake  reporting  capabilities  in  Fennoscandia  inferred  from  IAS  data  is 
given  in  Section  7.5,  while  Section  7.6  discusses  application  of  the  threshold 
monitoring  method  introduced  by  Ringdal  and  Kvaerna. 


2  NORSAR  Operation 


2.1  Detection  Processor  (DP)  Operation 

There  have  been  103  breaks  in  the  otherwise  continuous  operation  of  the  NOR¬ 
SAR  online  system  within  the  current  6-month  reporting  interval.  The  uptime 
percentage  for  the  period  is  91.8  as  compared  to  94.9  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  down¬ 
time  for  the  days  between  1  October  1989  and  31  March  1990.  The  monthly 
recording  times  and  percentages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  39 

b)  Stops  related  to  program  work  or  error  2 

c)  Hardware  maintenance  stops  18 

d)  Power  jumps  and  breaks  0 

e)  TOD  error  correction  31 

f)  Communication  lines  14 


The  total  downtime  for  the  period  was  338  hours  and  33  minutes.  The 
mean-time-between-failures  (MTBF)  was  1.6  days,  as  compared  to  1.9  for  the 
previous  period. 


J.  Torstveit 
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Fig.  2.1.1  Detection  Processor  downtime  for  October  (top),  November  (mid¬ 
dle)  and  December  (bottom)  1989. 


Fig.  2.1.1  Detection  Processor  downtime  for  January  (top),  February  (mid¬ 
dle)  and  March  (bottom)  1990. 
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21 

11 

290 

12 

20 

13 

290 

15 

0 

15 

290 

15 

0 

15 

291 

0 

0 

6 

291 

7 

25 

7 

291 

8 

34 

9 

291 

10 

12 

12 

291 

12 

46 

13 

291 

17 

48 

24 

292 

0 

0 

6 

292 

7 

27 

24 

293 

0 

0 

12 

processing  the  last 

COMMENTS . 

10  HODCOHP  FAILURE 
9  MODCOKP  FAILURE 
50  MODCOMP  FAILURE 
58  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 

10  MODCOMP  FAILURE 

46  MODCOMP  FAILURE 
21  MODCOMP  FAILURE 

47  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 
5  MODCOMP  FAILURE 

45  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 

40  MODCOMP  FAILURE 
55  MODCOMP  FAILURE 
45  MODCOMP  FAILURE 

41  MODCOMP  FAILURE 
40  MODCOMP  FAILURE 

11  MODCOMP  FAILURE 
47  MODCOMP  FAILURE 
52  CE  MAINTENANCE 

58  MODCOMP  FAILURE 
7  MODCOMP  FAILURE 

29  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 
16  MODCOMP  FAILURE 
32  CE  MAINTENANCE 
21  MODCOMP  FAILURE 
40  MODCOMP  FAILURE 
23  CE  MAINTENANCE 
44  CE  MAINTENANCE 

20  MODCOMP  FAILURE 
26  MODCOMP  FAILURE 
16  MODCOMP  FAILURE 

37  MODCOMP  FAILURE 

16  CE  MAINTENANCE 
35  MODCOMP  FAILURE 

18  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 

21  MODCOMP  FAILURE 
0  CE  MAINTENANCE 

52  CB  MAINTENANCE 

59  MODCOMP  FAILURE 

38  MODCOMP  FAILURB 
0  MODCOMP  FAILURE 

17  MODCOMP  FAILURE 

31  MODCOMP  FAILURB 

42  CB  MAINTENANCE 

43  CB  MAINTENANCE 
43  CE  MAINTENANCE 
25  MODCOMP  FAILURB 

32  MODCOMP  FAILURB 

19  MODCOMP  FAILURB 
0  MODCOMP  FAILURE 
4  MODCOMP  FAILURE 
0  MODCOMP  FAILURE 

20  MODCOMP  FAILURE 
0  CE  MAINTENANCE 

58  CE  MAINTENANCE 


HALF-YEAR 

DAY 

299 

304 

313 

315 

315 

319 

325 

326 

326 

327 
332 

332 

333 
335 
335 
339 
339 
347 
347 
347 
355 
355 
362 

1 

1 

2 

3 

4 

5 
9 

10 

18 

23 

29 

34 

35 

36 

36 

37 

38 

39 

40 
43 

43 

44 

46 

47 

58 

59 
59 
66 
72 
80 
80 
80 
82 
85 

ll 


START 

STOP 

COMMENTS . 

7 

1 

7 

3 

TOO  RETIRED  22MS 

7 

1 

7 

3 

TOD  RETARED  20MS 

7 

1 

7 

3 

TOD  RETARED  2 QMS 

13 

16 

14 

16 

CE  MAINTENANCE  MODC 

14 

45 

15 

42 

CE  MAINTENANCE  MODC 

7 

0 

7 

2 

TOD  RETARED  23MS 

8 

0 

0 

2 

TOD  RETARED  1 3MS 

7 

0 

7 

2 

TOD  RETARED  I2MS 

10 

28 

10 

30 

CE  MAINTENANCE  02B 

11 

51 

11 

53 

LINE  FAILURE 

7 

0 

7 

2 

TOD  RETARED  I4MS 

14 

44 

14 

46 

LINE  PAILURE 

11 

47 

11 

49 

LINE  FAILURE 

7 

53 

7 

55 

LINE  FAILURE 

8 

54 

8 

57 

LINE  FAILURE 

7 

40 

7 

42 

TOD  RETARED  22MS 

7 

42 

9 

31 

HARDWARE  FAILURE 

7 

41 

7 

4  3 

TOD  RETARED  22MS 

14 

33 

14 

35 

LINE  FAILURE 

12 

2 

12 

4 

LINE  FAILURE 

7 

4 

7 

6 

TOD  RETARED  23MS 

12 

47 

13 

7 

CE  MAINTENANCE 

10 

3 

10 

5 

LINE  FAILURE 

13 

35 

14 

0 

SYSTEM  WORK  (NEW  YEAR) 

14 

28 

14 

33 

SYSTEM  WORK  (NEW  YEAR) 

7 

25 

7 

27 

TOD  RETARED  90MS 

6 

30 

16 

30 

LINE  FAILURE 

10 

18 

18 

39 

LINE  FAILURE 

7 

1 

7 

3 

TOD  RETARED  13MS 

13 

10 

15 

45 

LINE  FAILURE 

14 

15 

14 

19 

CE  MAINTENANCE 

8 

51 

8 

53  LINE  FAILURE 

7 

I 

7 

3 

TOD  RETARED  2 2 MS 

7 

2 

7 

4 

TOD  RETARED  23MS 

8 

50 

24 

0 

LINE  FAILURE 

0 

0 

24 

0 

LINE  FAILURE 

0 

0 

10 

0 

LINE  FAILURE 

19 

40 

24 

0 

TOD  FAILURE 

0 

0 

24 

0 

TOD  FAILURE 

0 

0 

16 

48 

TOD  FAILURE 

12 

30 

12 

34 

TOD  RETARED  IONS 

7 

35 

9 

49 

MAINTENANCE  DISK 

6 

32 

6 

42 

TOD  FAILURE 

9 

31 

9 

43 

MAINTENANCE  TOD 

7 

10 

7 

12 

TOD  RETARED  12MS 

10 

6 

10 

8 

LINE  FAILURE 

7 

1 

7 

3 

TOD  RETARED  12MS 

7 

2 

7 

4 

TOD  RETARED  23MS 

8 

18 

8 

25 

TOD  FAILURE 

14 

45 

14 

48 

TOD  FAILURE 

7 

1 

7 

3 

TOD  RETARED  23MS 

7 

2 

7 

5 

TOD  RETARED  46MS 

7 

0 

7 

3 

TOD  RETARED  23MS 

10 

40 

10 

46 

TOD  FAILURE 

11 

3 

11 

6  TOD  FAILURE 

7 

0 

15 

10 

MODCOMP  FAILURE 

f 

11 

11 

36 

MODCOMP  FAILURE 

i 

44 

10 

8 

TOD  FAILURE 

10 

55 

11 

5 

TOD  FAILURE 

7 

47 

7 

49 

TOD  RETARED  3 MS 

Table  2.1.1  Daily  DP  downtime  in  the  period  1  October  1989  -  31  March 
1990. 
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Month  DP  Uptime  DP  Uptime  No.  of  No.  of  Days  DP  MTBF* 


Hours 

% 

DP  Breaks 

with  Breaks 

(days) 

OCT  89 

493.53 

66.3 

19 

NOV  89 

717.75 

99.7 

11 

8 

2.5 

DEC  89 

741.57 

99.7 

10 

6 

3.4 

JAN  30 

722.20 

97.1 

11 

10 

2.5 

FEB  90 

574.73 

85.5 

12 

13 

1.8 

MAR  90 

729.63 

98.1 

10 

7 

2.8 

91.1 

104 

63 

1.6 

♦Mean-time-between-failures  =  total  uptime/no.  of  up  intervals. 


Table  2.1.2  Online  system  performance,  1  October  1989  -  31  March  1990. 


2.2  Array  communications 
General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  through¬ 
out  the  reporting  period. 

TM-data  continued  the  fault  finding  procedure  on  the  Modcomp  computer 
in  October,  also  the  CPU  microcodes  were  analyzed. 

A  new  Modcomp  rack  equipped  with  new  boards  (except  for  2  memory 
boards  and  a  number  of  CPU  PROM  chips)  was  installed.  Also  this  initiative 
failed.  New  PROMs  were  programmed  and  installed  (23  Oct).  Since  that  date 
the  system  has  been  reliable. 

We  have  in  the  period  also  had  a  large  number  of  breaks  in  the  data 
stream  between  the  subarrays  and  NDPC.  Not  only  single  subarrays  have  been 
affected,  but  also  groups  of/ or  all  subarrays  have  been  affected  simultaneously. 

02C  was  affected  in  October  (week  40),  probably  caused  by  NTA  trouble. 
06C  outages  occurred  four  times  in  October  (weeks  40,  42,  43  and  44)  due  to 
stuck  CTV  equipment. 

In  November  a  new  break  in  the  data  stream  from  02B  occurred  (21  Nov). 
The  SLEM  Digital  Unit  (DCU)  was  replaced.  Also  02C  was  affected  (28  Nov) 
in  connection  with  scheduled  NTA  cable  work.  The  fined  operation,  splicing 
of  the  communications  cable,  went  wrong  and  resulted  in  a  break  towards  the 
CTV  for  3  days. 

7  December  was  the  start  of  a  bad  period  as  far  as  communications  was 
concerned.  NTA/Lillestram,  in  cooperation  with  NTA/Hamar  changes  to  an¬ 
other  carrier  system,  in  order  to  increase  the  capacity  between  Lillestrpm  and 
Hamar.  After  this  change  the  NORSAR  communications  system  became  less 
reliable.  All  systems  were  involved,  separately  and  simultaneously,  but  the 
systems  also  operated  well  over  several  days. 

The  communications  system  instability  continued  in  January.  Due  to  the 
special  character  of  the  error  pattern,  it  was  difficult  to  locate  the  source  for 
the  trouble.  Individual  end-to-end  subarray  communications  tests  during  bad 
periods  did  not  reveal  errors.  A  coax  carrier  system  common  to  all  systems 
functioned  well  until  7  January,  when  a  damaged  coax  carrier  cable  was  located 
in  Nittedal,  affecting  all  subarrays.  Apart  from  spikes  on  all  systems  12  Jan 
and  on  02C  and  03C  10  Jan,  the  system’s  performance  for  the  second  part  of 
January  improved. 

In  spite  of  2  damaged  coax  carrier  cables  located/repaired,  we  still  had 
irregularities  in  February.  The  two  last  outages  were  caused  by  the  Time 
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of  Day  Generator,  supplying  2400  Hz  to  the  NORSAR  systems.  Two  cards 
were  replaced  in  the  TOD,  and  the  situation  changed  radically  with  respect 
to  overall  system  performance,  although  individual  systems  were  affected  by 
other  problems  (02C,  03C  and  06C). 

In  March  the  Modcomp’s  disk  system  failed  two  times:  23  March  when 
the  Spindle  Drive  Motor  failed,  and  24  March  when  a  disk  “crash”  occurred, 
causing  a  total  outage  of  13  hrs  and  48  min. 

Also  the  TOD,  which  has  as  its  main  task  “time  tagging”  NORSAR  data 
and  clocking  out  the  NORSAR  commands  (ICWs),  has  been  less  reliable  than 
before.  The  oscillator  is  drifting  and  the  2400  Hz  has  been  “jittery”.  We  have 
made  an  appointment  with  a  firm  capable  of  handling  delicate  instruments 
such  as  the  TOD. 

The  performance  of  the  communications  systems,  including  communica¬ 
tions  lines,  carrier  systems,  modems,  etc.,  was  very  satisfactory  throughout 
March. 

Detailed  summary 

October  (weeks  40-44,  -41),  2.10-5.11.89 

Fault  finding  continued  with  the  Modcomp  by  TM-data.  CPU  microcodes 
were  analyzed,  and  irregularities  were  observed.  A  new  Modcomp  rack  equipped 
with  new  boards  (except  for  2  memory  boards  and  a  certain  number  of  PROM 
chips).  Even  then  the  system  failed,  and  new  PROMs  were  programmed 
and  installed  23  October.  The  Modcomp  system  has  been  reliable  after  that 
change. 

02C  was  down  week  40  (5.35%),  probably  due  to  NTA  trouble.  06C  out¬ 
ages  occurred  week  40  (7.14%),  week  42  (41.0%)  and  week  44  (12.5%)  due  to 
“stuck”  CTV  equipment. 

November  (weeks  45-48),  6.11-5.12.89 

The  Modcomp  was  reliable  after  repair  and  restart  23  October.  02B  was 
affected  21  November  due  to  a  failing  Digital  Unit  (DCU)  in  the  CTV.  23 
November  the  unit  was  replaced.  Scheduled  NTA  cable  work  north  of  Lille- 
hammer  affected  02C  28  November.  The  operation  involved  cable  cut  and 
splicing  afterwards.  The  latter  operation  failed,  causing  a  3-day  outage  to¬ 
wards  the  CTV. 

December  (weeks  49-52),  4-12-31.12.89 

7  December  NTA/Lillestrem  and  NTA/Hamar  in  cooperation  changed  a 
certain  carrier  group  in  order  to  increase  the  capacity  between  Lillestrom  and 
Hamar.  After  this  change  the  NORSAR  communications  system  became  less 


9 


reliable,  and  this  situation  continued  through  31  December  1989.  All  the  sub- 
array  communications  systems  were  affected,  both  simultaneously  and  sepa¬ 
rately.  The  systems  also  operated  satisfactorily  over  several  days. 

January  ( weeks  1~4),  1-28.1.90 

The  instability  which  started  in  December  continued  in  January.  It  was 
difficult  to  locate/isolate  the  source  of  the  trouble  due  to  the  special  character¬ 
istics  of  the  error  pattern.  Individual  end-to-end  tests  during  bad  periods  did 
not  reveal  irregularities/errors  (01A,  02C,  03C  and  02B  between  Lillestrom 
and  Hamar).  However,  symptoms  could  indicate  a  faulty  coax  carrier  ca¬ 
ble/system  common  to  all  NORSAR  systems. 

4  January  a  damaged  coax  cable  was  located  in  Lillestrom  in  connection 
with  construction  work.  After  the  repair  we  believed  the  problems  were  solved 
because  everything  functioned  properly  until  7  January.  This  time  a  damaged 
coax  carrier  cable  was  located  in  the  Nittedal  community,  affecting  all  the  sub¬ 
arrays.  2.5  hours  later  the  channels  carried  by  the  coax  had  been  transferred 
to  a  spare  coax  cable. 

Apart  from  spikes  on  all  channels  12  January,  and  02C/03C  10  January, 
the  system  performance  for  the  second  part  of  January  improved. 

February  f weeks  5-8),  29. 1-25. 2. 90 

Communications  system  instability  continued  in  February,  and  the  differ¬ 
ent  sources  of  trouble  these  weeks  are  as  follows: 

•  Damaged  coax  carrier  cable  Lillestrom 

•  Coax  carrier  cable  between  Kongsvinger  and  Skarnes  (reason  not  stated) 

•  Coax  carrier  cable  Nittedal  (gradual  deterioration) 

•  Time  of  Day  Generator  (TOD) 

The  two  last  outages  5-6  February  and  11-12  February  were  caused  by  the 
TOD;  two  cards  were  replaced. 

After  the  repair  the  situation  improved  radically  as  far  as  overall  system 
performance  is  concerned,  although  individual  systems  such  as  02C,  03C  and 
06C  have  been  affected  by  other  problems. 

March  (weeks  9-13),  26.2-1.4-90 

This  period  the  Modcomp  disk  system  failed  two  times.  23  March  the 
Spindle  Drive  Motor  stopped,  causing  an  outage  of  8  hrs  13  mins.  The  next 
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failure  occurred  24  March,  and  this  time  we  had  a  disk  “crash”,  causing  an 
outage  of  5  hrs  and  35  mins. 

Also  the  Time  Code  Generator  (TOD)  has  been  less  reliable.  The  oscillator 
is  drifting,  and  the  2400  Hz  clocking  the  NORSAR  commands  (IWUs)  toward 
the  CTV  has  been  “jittery”.  The  latter  may  have  been  caused  by  the  same 
“noise  spikes”  which  probably  also  causd  the  TOD  seconds  to  “jump”.  We 
believe  the  noise  spikes  observed  in  the  TOD  were  introduced  by  the  recorder 
used  during  time  checks.  For  these  checks  we  have  now  hooked  up  a  digital 
clock  locked  to  a  West  German  station. 

Otherwise  the  performance  of  the  communications  lines,  carrier  systems, 
modems  etc.,  have  been  most  satisfactory  throughout  the  month. 


O.A.  Hansen 


Sub¬ 

arrays 

Oct  89  (4) 
2-8.10  Sc 
16.10-5.11 

Nov  89  (4) 
6.11-3.12 

Dec  89  (4) 
4-31.12 

Jan  90  (4) 
1-28.1 

Feb  90  (4) 
29.1-25.2 

Mar  90  (5) 
26.2-1.4 

Average 
1/2  year 

01A 

*0.0040 

0.0030 

*0.500 

*1.058 

'“1*0.937 

0.013 

0.419 

01B 

0.0004 

0.0004 

*0.430 

*1.008 

">*0.936 

0.010 

0.397 

02B 

0.0020 

0.0340 

*0.440 

*1.059 

“>*0.926 

0.008 

0.411 

02C 

^•o.oooo 

3>*0.0006 

*0.447 

*0.580 

">*1.805 

0.009 

0.478 

03C 

0.0030 

0.0020 

*0.468 

*0.815 

">*2.696 

0.008 

0.615 

04C 

0.0030 

0.0006 

*0.462 

*0.974 

”>*0.900 

0.008 

0.391 

06C 

2>*  N/A 

0.0008 

*>*0.307 

5>*0.160 

">*1.187 

0.012 

7>*0.333 

AVER 

0.0020 

0.0060 

0.440 

0.807 

1.341 

0.009 

0.442 

LESS 

06C 

•See  Section  2.2  regarding  figures  preceded  by  an  asterisk. 

Figures  representing  error  rate  (in  per  cent)  preceded  by  a  number  1),  2),  etc.,  are  related  to 
legend  below. 


l>  Average  3  weeks  42-44 

2>  Stuck  CTV  equipment 

3>  Average  3  weeks  45-47,  scheduled  cable  work 

4>  Average  3  weeks 

5>  Average  2  weeks 

a>  Average  2  weeks 

7>  Average  5  months 


Dable  2.2.1  Communications  performance.  The  numbers  represent  error  rates 
in  per  cent  based  on  total  transmitted  frames/week  (1  October  1989  -  31 
March  1990). 
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2.3  Event  Detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor 
operation  are  given.  The  table  lists  the  total  number  of  detections  (DPX) 
triggered  by  the  on-line  detector,  the  total  number  of  detections  processed  by 
the  automatic  event  processor  (EPX)  and  the  total  number  of  events  accepted 
after  analyst  review  (teleseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 
P-phases  Core  phases 

Sum 

Daily 

Oct  89 

8525 

999 

183 

47 

230 

7.4 

Nov  89 

12575 

1574 

278 

53 

331 

11.0 

Dec  89 

18948 

2610 

260 

48 

308 

9.9 

Jan  90 

14624 

1622 

196 

64 

260 

8.4 

Feb  90 

10050 

1040 

199 

35 

234 

8.4 

Mar  90 

11750 

1209 

231 

54 

285 

9.2 

1347 

301 

1648 

9.1 

Table  2.3.1.  Detection  and  Event  Processor  statistics,  1  October  1989  -  31 
March  1990. 


B.  Paulsen 


3  Operation  of  regional  arrays 

3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 

Table  3.1.1  lists  the  main  outage  times  and  reasons,  and  as  can  be  seen  the 
main  reason  for  the  outage  is  hardware  failure  at  the  HUB  and  line  failure. 

The  average  recording  time  was  96.9%  as  compared  to  92.1%  for  the  pre¬ 
vious  period. 


Date 

Time 

Cause 

2 

Oct 

1756-1844 

Transmission  line 

failure 

7 

Oct 

0459-0906 

Transmission  line 

failure 

17 

Oct 

2306- 

Transmission  line 

failure 

18 

Oct 

-1253 

Transmission  line 

failure 

25 

Oct 

0336-0408 

Transmission  line 

failure 

25 

Oct 

0440-0628 

Transmission  line 

failure 

2 

Hov 

2346-2358 

Transmission  line 

failure 

12 

Hov 

2248- 

Hardware  failure  HUB 

13 

Hov 

- 

Hardware  failure  HUB 

14 

Hov 

-1149 

Hardware  failure  HUB 

22 

Dec 

1447- 

Hardware  failure  HUB 

23 

Dec 

-0806 

Hardware  failure  HUB 

24 

Dec 

0201-0603 

Hardware  failure  HUB 

28 

Dec 

0658-0808 

Transmission  line 

failure 

3 

Jan 

1206-1227 

Transmission  line 

failure 

4 

Jan 

1018-1843 

Transmission  line 

failure 

5 

Jan 

1625-1904 

Transmission  line 

failure 

6 

Jan 

2143- 

Transmission  line 

failure 

6 

Jan 

-1019 

Transmission  line 

failure 

11 

Jan 

2033-2108 

Transmission  line 

failure 

12 

Jan 

0944-1119 

Transmission  line 

failure 

16 

Jan 

0612-0617 

Transmission  line 

failure 

18 

Jan 

0943-0954 

Hardware  failure  HDPC 

25 

Jan 

0758-0927 

Transmission  line 

test 

31 

Jan 

1206-1212 

Hardware  failure  HDPC 

17 

Feb 

0114-0117 

Transmission  line 

failure 

18 

Feb 

0909-0912 

Transmission  line 

failure 

7 

Mar 

1300- 

Transmission  line 

test 

8 

Mar 

-1335 

Transmission  line 

test 

14 


13 

Mar 

1046-1108 

13 

Mar 

1606-1802 

13 

Mar 

1966-2106 

26 

Mar 

0712  0734 

26 

Mar 

1760  1834 

31 

Mar 

0640-0734 

31 

Mar 

2237-2334 

Hardware  failure  HUB 
Hardware  failure  HUB 
Hardware  failure  HUB 
Hardware  failure  HUB 
Hardware  failure  HUB 
Power  break  HUB 
Power  break  HUB 


Table  3.1.1.  Interruptions  in  NORESS  recordings  at  NDPC,  October  1989 
-  March  1990. 

Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 


October 

97.1% 

November 

94.8% 

December 

97.7% 

January 

96.2% 

February 

99.9% 

March 

95.9% 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 
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Fig.  3.1.1.  NORESS  data  recording  uptime  for  October  (top),  November 
(middle)  and  December  (bottom)  1989. 
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Fig.  3.1.1.  NORESS  data  recording  uptime  for  January  (top),  February 
(middle)  and  March  (bottom)  1990. 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjeller 


The  main  reasons  causing  most  of  the  ARCESS  outage  in  the  period  are: 
Hardware  failure  at  NDPC  or  HUB,  transmission  failure  and  power  failure  at 
HUB. 

The  long  outage  in  March  was  caused  by  clock  problems  at  HUB. 

The  average  recording  time  was  95.3%  as  compared  to  92.2%  for  the  pre¬ 
vious  period. 


Date 

Tine 

Cause 

10  Oct 

0242-0658 

Hardware  failure  at  NDPC 

12  Oct 

0813-0828 

Hardware  maintenance  at  NDPC 

13  Oct 

0954-1002 

Transmission  line  failure 

3  Hot 

1145-1168 

Transmission  line  failure 

20  Hot 

0504-0549 

Software  failure  at  NDPC 

28  Hot 

1110-1124 

Transmission  line  failure 

9  Dec 

0738-0825 

Hardware  failure  at  NDPC 

24  Dec 

0431- 

Powerbreak  HUB 

25  Dec 

- 

Powerbreak  HUB 

26  Dec 

-1027 

Powerbreak  HUB 

29  Dec 

0134-1707 

Hardware  failure  at  NDPC 

29  Dec 

2136- 

Hardware  failure  at  NDPC 

30  Dec 

-0610 

Hardware  failure  at  NDPC 

4  Jan 

2004-2218 

Hardware  failure  at  NDPC 

9  Jan 

0660-0946 

Transmission  line  failure 

12  Jan 

2044-2067 

Hardware  failure  at  NDPC 

14  Jan 

0938-1006 

Transmission  line  failure 

23  Jan 

1421-1443 

Transmission  line  failure 

28  Jan 

1700-1800 

Transmission  line  failure 

9  Feb 

0349-0407 

Hardware  failure  at  NDPC 

2306- 


15  Mar 

16  Mar 

17  Mar 

18  Mar 

19  Mar 


-2200 


21  Mar  1431-1906 
25  Mar  0800-1319 


27  Mar  1630-1707 


Hardware  failure  at  HUB 
Hardware  failure  at  HUB 
Hardware  failure  at  HUB 
Hardware  failure  at  HUB 
Hardware  failure  at  HUB 
Powerbreak  HUB 
Powerbreak  HUB 
Hardware  failure  at  NDPC 


Table  3.2.1  The  main  interruptions  in  ARCESS  recordings  at  NDPC,  Octo¬ 
ber  1989  -  March  1990. 

Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 

October  :  99.3% 

November  :  99.8% 

December  :  89.4% 

January  :  99.0% 

February  :  99.9% 

March  :  84.5% 

Fig.  3.2.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  October  (top),  November 
(middle)  and  December  (bottom)  1989. 


Fig.  3.2.1.  ARCESS  data  recording  uptime  for  January  (top),  February 
(middle)  and  March  (bottom)  1990. 
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3.3  Recording  of  FINESA  data  at  NDPC,  Kjeller 

Regular  recording  at  NDPC,  Kjeller,  of  data  from  the  FINESA  regional  array 
in  Finland  started  on  1  January  1990.  A  description  of  the  FINESA  system 
was  given  in  our  Semiannual  Technical  Summary  for  the  period  1  April  -  30 
September  1989.  The  average  recording  time  was  98.5%. 


Date 

Time 

Cause 

2  Jan 

0622-0722 

Hardware  failure  at  NDPC 

2  Jan 

0763-0823 

Transmission  line  failure 

11  Jan 

2033-2108 

Transmission  line  failure 

IS  Jan 

0759- 

Transmission  line  failure 

16  Jan 

-1261 

Transmission  line  failure 

1  Fab 

0002-0033 

Transmission  line  failure 

Table  3,3.1  The  main  interruptions  in  FINESA  recordings  at  NDPC,  January 
-  March  1990. 

Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 

January  :  95.8% 

February  :  99.9% 

March  :  99.9% 


Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  FINESA  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 
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Fig.  3.3.1.  FINESA  data  recording  uptime  for  January  (top),  February  (mid¬ 
dle)  and  March  (bottom)  1990. 
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3.4  Event  detection  operation 
NORESS  detections 

The  number  of  detections  (phases)  reported  during  day  274,  1989  through 
day  090,  1990  was  44113,  giving  an  average  of  244  detections  per  processed  day 
(181  days  processed). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  274,  1989  through  090,  1990,  2076  local  and  regional  events  were 
located  by  NORESS,  based  on  automatic  association  of  P-  and  S-type  arrivals. 
This  gives  an  average  of  11.5  events  per  processed  day  (181  days  processed).  67% 
of  these  events  are  within  300  km,  and  88%  of  these  events  are  within  1000  km. 

ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  274,  1989  through 
day  090,  1990  was  66702,  giving  an  average  of  375  detections  per  processed  day 
(178  days  processed). 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  274, 1989,  through  090, 1990, 3144  local  and  regional  events  were 
located  by  ARCESS,  based  on  automatic  association  of  P-  and  S-type  arrivals. 
This  gives  an  average  of  16.9  events  per  processed  day  (178  days  processed).  54% 
of  these  events  are  within  300  km,  and  87%  of  these  events  are  within  1000  km. 

FINESA  detections 

The  number  of  detections  (phases)  reported  during  day  327,  1989  through 
day  090,  1990  was  29240,  giving  an  average  of  290  detections  per  processed  day 
(101  days  processed). 

Table  3.4.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 

Events  automatically  located  by  FINESA 

During  days  327, 1989,  through  090, 1990,  2000  local  and  regional  events  were 
located  by  FINESA,  based  on  automatic  association  of  P-  and  S-type  arrivals. 
This  gives  an  average  of  20.4  events  per  processed  day  (98  days  processed).  74% 
of  these  events  are  within  300  km,  and  89%  of  these  events  are  within  1000  km. 

J.  Fyen 
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each  hour  of  the  day  and  number  of  detections  for  that  day.  The  end  statistics  give  total  numDer  01  detections  aisinouiea 
for  each  hour  and  the  total  sum  of  detections  during  the  period.  The  averages  show  number  of  processed  days,  hourly 
distribution  and  average  per  processed  day.  (Page  4  of  4) 
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Table  3.4.3  (Page  1  of  3) 


3.5  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  De¬ 
cember  1989  and  has  been  operated  experimentally  since  1  January  1990  for 
automatic  processing  of  multiple-array  data.  The  current  version  of  IMS  pro¬ 
cesses  data  from  the  two-array  network  consisting  of  NORESS  and  ARCESS. 
Future  upgrades  of  IMS  will  allow  data  from  additional  arrays  and  single  sta¬ 
tions  to  be  incorporated. 

Table  3.5.1  gives  a  summary  of  phase  detections  and  processed  regional 
events  by  IMS  during  its  first  three  months  of  operation  at  NORSAR.  From 
top  to  bottom,  the  table  gives  the  total  number  of  detections  by  the  IMS,  the 
detections  that  are  associated  with  regional  events  declared  by  the  IMS,  the 
number  of  detections  that  are  not  associated  with  such  events,  the  number  of 
regional  events  declared  by  the  IMS,  the  number  of  such  events  rejected  by  the 
analyst,  the  number  of  events  accepted  by  the  analyst,  the  number  of  events 
accepted  by  the  analyst  without  any  changes,  and  finally  the  number  of  events 
accepted  after  some  sort  of  modification  by  the  analyst.  This  last  category 
is  divided  into  three  classes:  Events  where  phases  (not  detected  by  the  IMS) 
have  been  added  by  the  analyst,  events  for  which  the  phase  assignments  by  the 
IMS  have  been  changed  or  one  or  more  phase  detections  have  been  removed, 
and  events  for  which  the  changes  by  the  analyst  have  amounted  to  retiming 
the  phases  only. 


B.Kr.  Hokland 
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Jan  90 

Feb  90 

Mar  90 

total 

Phase  detections 

15421 

15313 

15012 

45746 

-  Associated  phases 

2043 

2485 

2244 

6772 

-  Unassociated  phases 

13378 

12828 

12768 

38974 

Events  declared 

1131 

1174 

982 

3287 

-  Rejected  events 

403 

355 

234 

992 

-  Accepted  events 

728 

819 

748 

2295 

Unchanged  events 

261 

127 

124 

512 

Modified  events 

467 

692 

624 

1783 

Phases  added 

39 

49 

64 

152 

Phases  changed  or  removed 

318 

437 

379 

1134 

Retiming  only 

110 

206 

181 

497 

Table  3.5.1.  IMS  phase  detections  and  events  summary. 


4  Improvements  and  modifications 

4.1  NORSAR 

NORSAR  Data  acquisition 

No  modification  has  been  done  to  the  NORSAR  data  acquisition  system. 

The  process  of  evaluating  technical  options  for  upgrading  the  array  is  con¬ 
tinuing.  The  requirement  is  to  find  digitizer  units  that  will  fit  into  existing 
seismometer  vaults,  and  that  may  use  existing  buried  cables  for  DC  power 
and  data  transmission  to  each  subarray  center.  So  far  it  seems  to  be  possible 
to  use  the  Nanometric  RD3  digitizers  at  the  seismometer  with  help  of  small 
batteries.  These  digitizers  give  24-bit  dynamic  range  using  16-bit  resolution 
and  automatic  gain  ranging  (comparable  to  NORESS).  Digitizers  with  24-bit 
resolution  are  also  available,  but  the  power  requirement  for  most  existing  units 
is  currently  beyond  the  supported  level  at  the  seismometer  site. 

Continued  research  will  be  done  to  find  24-bit  resolution  digitizers  with 
satisfactory  power  consumption. 

Data  center  operation,  communication  with  subarrays,  and  subarray  pro¬ 
cessing  systems  in  connection  with  the  upgrade  appear  to  pose  few  problems, 
as  many  of  the  NORESS/FINESA/GERESS/Hagfors  concepts  may  be  used. 

A  test  of  a  full  subarray  acquisition  system  will  be  performed  at  the  end 
of  the  next  reporting  period. 

NORSAR  Detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the 
IBM  4341/4381  computers  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2. 

NORSAR  Event  processing 

There  are  no  changes  in  the  routine  processing  of  NORSAR  events,  using 
the  IBM  system.  Full  event  processing  on  the  SUN  system  is  not  yet  imple¬ 
mented.  Processing  of  the  large  NORSAR  array  requires  a  data  base  of  time 
delay  corrections  to  get  the  benefit  of  using  the  array  for  locating  teleseismic 
events.  The  current  time  delay  correction  data  base  on  IBM  is  not  easy  to 
convert  to  SUN/unix  installations.  In  light  of  a  future  upgrade  of  the  system, 
a  research  task  would  be  to  build  a  new  time  delay  correction  data  base  for 
the  NORSAR  array,  with  a  design  that  may  be  used  for  any  type  of  array  and 
computer  system. 
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4.2  Regional  arrays 

Detection  processing 

The  routine  detection  processing  of  the  regional  arrays  is  running  satis¬ 
factorily  on  each  of  the  arrays’  SUN-3/280  acquisition  systems.  The  same 
program  is  used  for  NORESS/ARCESS/  FINESA,  but  with  different  recipes. 
The  beam  table  for  NORESS/ARCESS  may  be  found  in  NORSAR  Sci.  Rep. 
No.  1-89/90.  The  beam  deployment  for  FINESA  is  given  in  Table  4.2.1. 

Event  processing.  Phase  estimation 

This  process  applies  F-K  and  polarization  analysis  to  each  detection  in 
order  to  determine  phase  velocity,  azimuth  and  type  of  phase. 

The  processing  makes  use  of  the  EP  program,  and  may  be  performed  on 
any  of  the  SUN  workstations. 

Plot  and  epicenter  determination 

A  description  of  single  array  event  processing  has  been  given  in  NORSAR 
Sci.  Rep.  No.  2-88/89. 

The  processing  makes  use  of  the  EP  program,  and  may  be  performed  on 
any  of  the  SUN  workstations.  Different  recipes  are  used  for  each  array,  as  the 
rules  for  associating  phases  are  somewhat  different.  In  particular,  the  relation 
between  Sn,  Lg  and  Rg  in  terms  of  amplitudes  is  quite  different  for  the  three 
arrays. 

The  main  output  is  event  parameters  and  plots  of  single  sensor  and  beam 
traces  for  each  regional  and  teleseismic  event.  These  plots  are  useful  both  for 
assessing  processing  performance  and  in  monitoring  array  data  quality. 

The  routine  processing  of  FINESA  data  at  NORSAR  is  similar  to  what  is 
done  in  Helsinki. 

Automatic  single  array  event  bulletins  are  available,  but  the  solutions  are 
not  routinely  reviewed,  as  this  is  the  task  for  the  IMS  system. 


4.3  GSETT-2  experiments 

EP  recipes  for  extracting  NORESS/ARCESS  events  and  producing  level  1  and 
level  2  data  files  have  been  developed. 


The  EP  program  read*  phase  information  from  the  automatic  (or  reviewed) 
bulletins,  accesses  the  necessary  seismometer  data,  forms  array  beams,  and 
produces  level  2  data  files, 
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BEAN  Velocity  Azimuth  Filter  band  Threshold  Configuration 


coil 

99999.9 

0.0 

0.50 

- 

1.50 

4.20 

7 

A1 

C 

C021 

99999.9 

0.0 

1.0 

- 

3.0 

4.20 

7 

A1 

C 

C031 

99999.9 

0.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C032 

11.0 

30.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C033 

11.0 

90.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C034 

11.0 

150.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C03S 

11.0 

210.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C036 

11.0 

270.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C037 

11.0 

330.0 

1.50 

- 

3.50 

4.20 

7 

A1 

c 

C038 

15.0 

80.0 

1.50 

- 

3.50 

3.70 

7 

A1 

c 

C039 

10.0 

20.0 

1.50 

- 

3.50 

3.70 

7 

A1 

c 

C041 

99999.9 

0.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C042 

10.20 

30.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C043 

10.20 

90.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C044 

10.20 

150.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C045 

10.20 

210.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C046 

10.20 

270.0 

2.0 

- 

4.0 

4.20 

7 

Ai 

c 

C047 

10.20 

330.0 

2.0 

- 

4.0 

4.20 

7 

A1 

c 

C048 

15.0 

80.0 

2.0 

- 

4.0 

3.70 

7 

Al 

c 

C049 

10.0 

20.0 

2.0 

- 

4.0 

3.70 

7 

Al 

c 

C0S1 

99999.9 

0.0 

2.50 

- 

4.50 

4.20 

12 

Bi 

BC 

C052 

8.90 

30.0 

2.50 

- 

4.50 

4.20 

12 

B1 

BC 

C053 

8.90 

90.0 

2.50 

- 

4.50 

4.20 

12 

Bl 

BC 

C054 

8.90 

150.0 

2.50 

- 

4.50 

4.20 

12 

Bl 

BC 

C0S5 

8.90 

210.0 

2.50 

- 

4.50 

4.20 

12 

Bl 

BC 

C056 

8.90 

270.0 

2.50 

- 

4.50 

4.20 

12 

Bl 

BC 

C0S7 

8.90 

330.0 

2.50 

- 

4.50 

4.20 

12 

Bl 

BC 

C0S8 

15.0 

80.0 

2.50 

- 

4.50 

3.70 

12 

Bl 

BC 

C061 

99999.9 

0.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C062 

10.50 

30.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C063 

10.60 

90.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C064 

10.50 

150.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C066 

10.50 

210.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C066 

10.50 

270.0 

3.0 

- 

6.0 

4.20 

12 

Bl 

BC 

C067 

10.50 

330.0 

3.0 

- 

5.0 

4.20 

12 

Bl 

BC 

C068 

15.0 

80.0 

3.0 

- 

5.0 

3.70 

12 

Bl 

BC 

C071 

99999.9 

0.0 

3.50 

- 

5.50 

4.20 

IS 

Al 

ABC 

C072 

11.10 

30.0 

3.50 

- 

5.50 

4.20 

15 

Al 

ABC 

C073 

11.10 

90.0 

3.50 

- 

5.50 

4.20 

15 

Al 

ABC 

C074 

11.10 

150.0 

3.50 

- 

6.60 

4.20 

15 

Al 

ABC 

C075 

11.10 

210.0 

3.50 

- 

5.50 

4.20 

15 

Al 

ABC 
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C076 

11.10 

270.0 

3.50 

- 

5.50 

4.20 

15 

A1 

ABC 

C077 

11.10 

330.0 

3.50 

- 

5.50 

4.20 

15 

A1 

ABC 

C081 

99999.9 

0.0 

4.0 

- 

8.0 

4.20 

15 

A1 

ABC 

C082 

9.50 

30.0 

4.0 

- 

8.0 

4.20 

15 

A1 

ABC 

C083 

9.50 

90.0 

4.0 

- 

8.0 

4.20 

15 

A1 

ABC 

C084 

9.50 

150.0 

4.0 

- 

8.0 

4.20 

IS 

A1 

ABC 

C085 

9.50 

210.0 

4.0 

- 

8.0 

4.20 

15 

A1 

ABC 

C086 

9.50 

270.0 

4.0 

- 

8.0 

4.20 

15 

A1 

ABC 

C087 

9.50 

330.0 

4.0 

- 

8.0 

4.20 

IS 

A1 

ABC 

C091 

99999.9 

0.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C092 

10.50 

30.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C093 

10.50 

90.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C094 

10.50 

150.0 

5.0 

- 

10.0 

4.70 

IS 

A1 

ABC 

C095 

10.50 

210.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C096 

10.50 

270.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C097 

10.50 

330.0 

5.0 

- 

10.0 

4.70 

15 

A1 

ABC 

C101 

99999.9 

0.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

C102 

9.90 

30.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

C103 

9.90 

90.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

C104 

9.90 

150.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

C105 

9.90 

210.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

C106 

9.90 

270.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

0107 

9.90 

330.0 

8.0 

- 

16.0 

4.70 

9 

A1 

AB 

CI01 

99999.9 

0.0 

0.50 

- 

1.50 

2.50 

7 

A1 

C 

CI02 

99999.9 

0.0 

1.0 

- 

2.0 

2.50 

7 

A1 

C 

CI03 

99999.9 

0.0 

1.50 

- 

2.50 

2.50 

7 

A1 

C 

CI04 

99999.9 

0.0 

2.0 

- 

4.0 

2.10 

7 

A1 

C 

CI05 

99999.9 

0.0 

3.50 

- 

5.50 

2.10 

7 

A1 

C 

CI06 

99999.9 

0.0 

5.0 

- 

10.0 

2.50 

7 

A1 

C 

Table  4.2.1.  FINESA  beam  table.  The  table  shows  the  name  of  the  beam, 
velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz),  STA/LTA  threshold, 
and  configuration.  The  configuration  is  described  with  number  of  sensors  and 
A,B  or  C-rings,  with  center  A1  or  Bl.  CI01  -  CI06  are  incoherent  beams. 
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5  Maintenance  Activities 


5.1  Activities  in  the  held  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  in  the  field,  at  the  Maintenance  Center 
(NMC)  Hamar  and  NDPC  activities  related  to  monitoring  and  control  of  the 
NORSAR,  NORESS  and  ARCESS  arrays. 

Activities  comprise  preventive/corrective  maintenance  in  connection  with 
all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  In  addition,  NMC  has 
been  involved  in  modification  of  equipment  (FINESA)  and  preparatory  work  in 
connection  with  HF  instrumentation  (NORESS).  Other  activities  are  related 
to  testing  the  NORSAR  communications  systems. 

NORSAR 

In  the  NORSAR  array,  spike  complaints  have  been  investigated,  SP/LP 
channel  gain  adjusted,  LP  seismometer  parameters  corrected,  SLEMs  have 
been  replaced,  batteries  and  transmitters  replaced  at  the  telemetry  stations, 
and  communications  lines  have  been  tested,  including  modems  and  level  checks. 

NORESS 

The  NORESS  array  has  been  visited  in  connection  with  preparing  vault  C2 
for  HF  instrumentation  instead  of  AO.  Hub  69  cards  (I/F)  have  been  replaced, 
and  the  Hub  power  supply  repaired.  NORESS  has  been  visited  in  connection 
with  communication  trouble,  all  equipment  in  vault  C2  was  removed  as  a 
precaution,  ADC  card  replacements  were  done. 

ARCESS 

ARCESS  activities  comprise  resetting  of  the  “white  box”,  replacement 
of  GPS  synchronized  clock  (satellite  clock),  and  replacement  of  fiber  optic 
transmitters. 


Details  are  presented  in  Table  5.1. 


Subarray/ 

area 

Task 

Date 

NORSAR: 

01A 

SA  visited  in  connection  with  spike  complaints 

18  Oct 

06C 

All  SP  channel  gains  adjusted. 

18  Oct 

NORESS: 

Vault  C2  made  ready  in  order  prepare  this  vault 
for  HF  instrumentation  instead  of  AO 

4-6  Oct 

Hub  69  card  (I/F  card)  repaired  due  to  a  bad  line 
driver 

18  Oct 

FINESA: 

Field  maintenance 

week  48 

NMC: 

General  work,  incl.  the  FINESA  acquisition  system 
equipment 

Oct 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 
Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Oct 

NORSAR: 

06C 

Adjusted  SP/LP  instruments. 

9  Nov 

02B 

Visit  in  connection  with  degraded  data. 

22  Nov 

04C/06C 

Adjustment  of  SP/LP  instruments. 

23  Nov 

02B 

Replaced  Digital  Unit  (DCU) 

24  Nov 

01B/02C 

Adjustment  of  SP/LP  instruments. 

27  Nov 

03C 

Adjustment  of  SP/LP  instruments. 

30  Nov 

NORESS: 

Hub  reset. 

21  Nov 

FINESA: 

Installation  of  a  new  data  acquisition 
system. 

13-16  Nov 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 
Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Nov 

Table  5.1  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  in¬ 
cluding  NDPC  activities  related  to  the  NORSAR  array,  1  October  1989  -  31 
March  1990. 


Subarray/ 

area 

Task 

Date 

NORSAR: 

02B(tel) 

Replaced  batteries  at  stations  1,  5  and  6. 

Replaced  transmitter  at  1  and  5. 

1,4, 7, 8 

9,12  and 

18  Dec 

01A 

Line  test  and  SP/LP  adjustment. 

6  Dec 

NORESS: 

ARCESS: 

Hub  power  supply  repaired. 

NORSAR’s  local  contact  visited  the  site  and 
reset  the  “white  box” . 

24  Dec 

FINESA: 

Data  acquisition  unit  returned  from  Finland  and 
modified  with  respect  to  hardware  and  software 

14  Dec 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 

Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Dec 

NORSAR: 

01A.06C 

Communications  test. 

3  Jan  1990 

NORESS: 

Visit  in  connection  with  communications  trouble 

6  Jan 

FINESA: 

Faulty  cable  located  between  B2  and  C5. 

Replaced  amplifier  (C6). 

30-31  Jan 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 

Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Jan 

NORSAR: 

01A 

Line  and  communications  tests. 

7  Feb 

06C 

Line  and  communications  tests. 

7  Feb 

06C 

NTA/Hamar  was  advised  to  check  the  line  level 
toward  the  CTV.  As  a  result,  the  level  was  raised 
by  5  dBm. 

15  Feb 

02B(tel) 

A  new  battery  was  installed  at  the  receiving  station. 

9  Feb 

Table  5.1  (cont.) 


Stibarray/ 

area 

Task 

Date 

NORESS: 

As  a  precaution,  all  equipment  in  the  vault  C2  was 
disconnected  and  removed. 

23  Feb 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 
Continuous  measurement  of  Mass  Position  and  Free 
Period 

Adjustment  of  Mass  Position  and  ftee  Period  when 
outside  tolerances 

Feb 

NORSAR: - - 

02B(tel) 

Replaced  station  4  battery. 

12  Mar 

ARCESS: 

Replaced  GPS  synchronized  clock  (satellite 
clock)  with  a  4F-DC  Radio  Clock. 

Replaced  fiber  optic  transmission  at  remote  site 

D6. 

Replaced  fiber  optic  transmission  at  Hub  site  B3. 

19-23  Mar 

NDPC: 

Daily  check  of  data  and  comm,  systems  NORSAR, 
NORESS,  ARCESS  and  FINESA. 

Weekly  calibration  of  SP/LP  instruments. 

Continuous  measurement  of  Mass  Position  and  fVee 
Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Mar 

Table  5.1  (cont.) 


5.2  Array  status 

As  of  31  Mar  1990  the  following  NORSAR  channels  deviated  from  tolerances. 

01A  01  8  Hz  filter 
02  8  Hz  filter 
04  30  dB  attenuator 
02C  04  Spikes 


O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Pub¬ 
lished 


7.1  Report  from  the  symposium  on  “Regional  Seismic  Ar¬ 
rays  and  Nuclear  Test  Ban  Verification”  in  Oslo,  Norway, 
14-17  February  1990 

During  14-17  February  1990  NORSAR  hosted  an  international  symposium 
entitled  “Regional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification”.  The 
symposium  was  attended  by  76  scientists  and  representatives  from  21  coun¬ 
tries,  including  a  large  number  of  seismologists  participating  in  the  work  of  the 
Conference  on  Disarmament’s  Group  of  Scientific  Experts  (GSE)  in  Geneva. 

The  purpose  of  the  symposium  was  to  assess  the  state-of-the-art  of  research 
on  regional  seismic  arrays  and  associated  topics.  In  particular,  the  symposium 
focused  upon  the  advanced  regional  arrays  NORESS  and  ARCESS  in  Norway 
and  their  associated  data  processing  facilities,  in  the  light  of  the  potential  of 
such  arrays  to  provide  a  much  improved  monitoring  capability  for  a  future 
comprehensive  nuclear  test  ban  treaty.  During  a  three-day  scientific  sympo¬ 
sium,  a  number  of  presentations  were  given  on  topics  relevant  to  this  issue. 
A  special  session  was  devoted  to  summarizing  the  experience  and  discussing 
further  plans  for  the  on-going  international  GSE  experiment  (GSETT-2). 

In  this  paper,  we  give  a  brief  review  of  some  of  the  results  presented  during 
the  scientific  symposium.  A  list  of  all  the  presentations  is  provided  in  an 
appendix,  and  the  numbers  in  brackets  refer  to  this  list.  The  majority  of  the 
papers  have  been  submitted  for  publication  to  the  Bulletin  of  the  Seismological 
Society  of  America,  and  will  be  published  in  a  Special  Issue  of  the  Bulletin, 
scheduled  to  appear  in  the  fall  of  1990. 

Development  of  regional  arrays 

Reviews  of  recent  developments  with  regard  to  regional  seismic  arrays  are 
presented  for  NORESS  and  ARCESS  in  Norway  [1  j,  GERESS  in  the  Federal 
Republic  of  Germany  [2]  and  FINESA  in  Finland  [3].  Paper  [1]  summarizes  the 
design  considerations  leading  to  the  establishment  of  the  first  regional  array, 
NORESS,  and  describes  how  the  success  of  this  new  array  concept  motivated 
the  deployment  of  additional  arrays  of  this  type.  The  paper  documents  the 
basic  signal  processing  techniques  used  in  real-time  data  analysis  for  regional 
arrays,  and  demonstrates  the  excellent  detection  performance  of  such  arrays 
at  regional  distances  (less  than  2000  km).  It  is  shown  that  NORESS  and 
ARCESS  are  capable  of  detecting  seismic  events  of  magnitude  2.5  with  90 
per  cent  probability,  if  these  events  occur  within  1000  km  distance,  whereas 
global  teleseismic  networks  have  much  higher  event  detection  thresholds.  The 
FINESA  array  is  also  documented  to  have  an  excellent  performance  [3],  and 
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together,  these  three  arrays  are  capable  of  locating  weak  seismic  events  in 
Fennoscandia  very  accurately  (typically  to  within  10-20  km).  The  GERESS 
array  currently  under  development  shows  many  of  the  same  excellent  features 
[2],  and  will  contribute  further  to  an  excellent  regional  coverage  of  large  parts 
of  Northern  Europe. 

Processing  of  data  from  a  network  of  regional  arrays 

Recent  technological  advances  have  allowed  very  sophisticated  processing  tech¬ 
niques  to  be  applied  in  detecting,  locating  and  identifying  seismic  events  using 
a  network  of  seismic  arrays  and  single  stations,  and  this  is  highlighted  by  the 
development  of  the  Intelligent  Monitoring  System  (IMS)  [4], [5].  Two  of  the 
goals  for  this  system  are  (1)  to  demonstrate  the  monitoring  performance  and 
capability  of  the  system  for  small  events  at  regional  distances  and  (2)  to  ex¬ 
plore  the  promise  of  an  expert-systems  approach  for  providing  improved  mon¬ 
itoring  performance  as  experience  accumulates.  The  first  operational  version, 
described  in  [4],  processes  data  from  NORESS  and  ARCESS,  whereas  later 
versions  will  be  expanded  to  networks  including  both  arrays  and  single  sta¬ 
tions.  The  IMS  is  ambitious  in  exploring  and  integrating  many  new  computer 
technologies,  and  the  validity  of  the  concept  is  documented  in  an  evaluation 
of  its  initial  operational  performance  [5]. 

Signal  analysis  methods 

A  number  of  presentations  addressed  methods  for  processing  seismic  signals 
recorded  by  arrays  as  well  as  three-component  stations.  It  was  demonstrated 
that  both  types  of  stations  can  provide  information  very  useful  in  phase  identi¬ 
fication,  azimuth  estimation  and  estimating  the  apparent  velocity  of  detected 
phases.  From  theoretical  considerations  as  well  as  from  experimental  com¬ 
parison  [12], [26], [14]  arrays  are  shown  to  be  superior  in  this  regard  at  low 
signal-to-noise  ratios,  although  the  precision  e.g.  of  azimuth  estimates  is  in¬ 
fluenced  by  a  number  of  factors,  including  phase  type,  frequency  of  the  signal 
and  systematic  bias  caused  by  earth  heterogeneities  [14], [18],  [26].  A  very 
promising  approach,  discussed  in  [17]  is  that  of  joint  analysis  of  3-component 
and  array  data. 

Signal  detection  methods  are  discussed  in  several  papers.  In  [11],  a  system 
for  on-line  detection  and  signal  analysis  is  presented  as  applied  to  a  Soviet  3- 
component  station  in  Kazakhstan.  In  [13],  a  detection  technique  is  described 
using  NORESS  array  and  3-component  data.  A  statistical  approach,  using 
adaptive  techniques,  to  detection  processing  and  estimation  is  presented  in  [7] 
for  array  data  and  [15]  for  3-component  data.  A  new  approach  to  obtain  precise 
relative  location  estimates  of  seismic  events,  using  high  frequency  recordings, 
is  presented  in  [25]. 
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Source  identification 

Traditionally,  seismic  discrimination  research  has  focused  on  distinguishing 
between  earthquakes  and  underground  nuclear  explosions.  Under  a  Compre¬ 
hensive  Test  Ban  Treaty,  emphasis  will  be  on  detecting  and  identifying  weak 
seismic  events,  and  a  third  category,  large  chemical  explosions  for  industrial 
purposes  (e.g.  mining  work)  will  become  important  to  consider.  In  [10],  a 
very  promising  method  is  applied  to  NORESS  data  to  discriminate  between 
earthquakes  and  ripple-fired  quarry  blasts  (mining  events  consisting  of  several 
explosions  closely  grouped  in  space  and  time).  Using  spectral  characteristics  of 
the  signals,  an  “automatic”  discriminant  is  proposed  computing  the  likelihood 
that  ripple  firing  occurred  in  each  given  case. 

In  [8],  a  novel  approach  making  use  of  artificial  neural  networks  is  used  to 
develop  a  classification  procedure  between  earthquakes  and  mining  explosions. 
Also  in  this  approach,  the  spectral  characteristics  of  the  signals  form  the  basis 
for  the  discriminants.  The  neural  network  appears  to  improve  in  particular 
the  classification  of  outliers  in  the  population,  and  reduce  the  number  of  un¬ 
certain  events.  Application  of  neural  networks  in  improving  seismic  processing 
performance  is  also  addressed  in  [9]. 

Of  considerable  interest  for  source  identification  is  also  the  method  pro¬ 
posed  in  [16],  applying  transfer  functions  to  transform  e.g.  between  recordings 
of  presumed  single  explosions  and  ripple-fired  explosions,  and  also  between 
recordings  at  different  NORESS  sensors  for  a  given  event.  This  gives  promise 
to  improve  the  coherence  of  seismic  phases  recorded  at  an  array,  with  ensuing 
implications  for  improved  source  parameter  estimation.  In  [6],  a  case-based 
reasoning  approach  to  event  identification  is  discussed,  and  a  waveform  enve¬ 
lope  matching  technique  is  applied  to  a  set  of  Western  Norway  earthquakes 
and  explosions. 

Detection  thresholds  and  in-country  networks 

While  regional  arrays  were  originally  designed  to  enhance  the  capabilities  for 
detecting  and  characterizing  weak  seismic  events  at  regional  distances,  they 
have  also  been  found  very  effective  in  the  teleseismic  distance  range.  As  an 
example,  published  yields  of  Soviet  underground  nuclear  explosions  at  Semi- 
palatinsk  have  been  used  to  evaluate  the  NORESS  detection  threshold,  in 
terms  of  explosive  yield  for  events  at  this  test  site  [21].  The  threshold  for  de¬ 
tection  at  NORESS  is  estimated  to  be  as  low  as  0.1  kt,  assuming  full  coupling 
and  normal  noise  conditions.  It  is  pointed  out  that  NORESS  has  particularly 
favorable  conditions  for  detecting  small  events  from  this  test  site,  and  that  the 
seismic  identification  threshold  necessarily  will  be  higher  than  the  detection 
threshold. 

Data  from  new  Global  Seismic  Network  stations  in  the  Soviet  Union,  in¬ 
stalled  as  a  cooperative  project  between  American  and  Soviet  scientists,  have 


been  applied  in  several  studies  to  address  problems  relevant  to  an  in-country 
monitoring  network.  Seismic  noise  levels  at  tbese  stations  are  analyzed  in 
[19],  and  found  to  be  higher  than  at  NORESS  in  the  band  1-20  Hz,  with  max¬ 
imum  difference  ranging  from  7  to  25  dB,  depending  on  the  station.  However, 
significant  noise  reduction  can  be  achieved  by  borehole  deployment. 

Using  data  from  stations  in  the  USSR,  the  frequency-dependent  attenua¬ 
tion  of  regional  seismic  phases  has  been  studied  in  [22].  Attenuation  charac¬ 
teristics  are  found  to  be  similar  to  those  observed  in  Scandinavia,  but  with  an 
absolute  Pn  amplitude  almost  a  factor  of  2  higher  in  eastern  Kazakhstan  for 
a  fixed  Lg  magnitude. 

Recordings  of  Semipalatinsk  nuclear  explosions  at  the  new  Global  Seismic 
Network  stations  in  the  Soviet  Union,  together  with  data  from  stations  in 
China  have  been  analyzed  in  [20]  and  it  is  shown  that  RMS  Lg  can  be  measured 
at  widely  separated  stations  with  a  remarkable  degree  of  consistency.  The 
standard  deviation  of  the  differences  between  pairs  of  stations  is  as  low  as 
0.03-0.04  in  logarithmic  units,  and  reliable  measurements  may  be  made  at 
magnitude  (mj)  down  to  about  4.0  for  stations  situated  about  1500  km  away 
from  Semipalatinsk.  The  importance  of  this  observation  in  terms  of  supplying 
yield  estimates  for  nuclear  explosions  down  to  and  even  below  one  kiloton  is 
pointed  out. 

Earth  structure,  wave  propagation,  scattering 

Several  of  the  papers  were  devoted  to  studies  of  general  problems  in  seismology 
and  geophysics,  in  areas  relevant  to  the  seismic  monitoring  issue.  The  struc¬ 
ture  of  the  crust  and  upper  mantle  in  parts  of  Northern  Eurasia  is  addressed 
in  papers  [23],  [24],  [27]  and  [29],  with  the  three  latter  papers  specifically  mak¬ 
ing  use  of  regional  array  data.  Seismic  wave  propagation  and  scattering  are 
addressed  in  a  number  of  papers,  e.g.  [13],  [26],  [28],  [29],  [30]. 

Conclusion 

The  Oslo  symposium  demonstrated  the  considerable  progress  in  the  field  of 
seismic  monitoring  during  recent  years.  It  particularly  highlighted  the  tech¬ 
nological  advances  in  seismic  instrumentation,  data  communication  and  com¬ 
puter  processing,  as  exemplified  by  the  development  of  advanced  regional  seis¬ 
mic  arrays  with  very  sophisticated  automatic  and  interactive  signal  process¬ 
ing  facilities.  The  presentations  at  the  scientific  symposium  show  that  these 
technological  advances  are  accompanied  by  considerable  scientific  progress,  al¬ 
though  much  work  remains  in  order  to  fully  exploit  the  potential  offered  by 
regional  arrays  in  a  seismic  monitoring  context. 


F.  Ringdal 
S.  Mykkeltveit 
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Appendix 


List  of  scientific  presentations  given  during  the  1990  Oslo  Symposium  on 
Regional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification 
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7.2  Analysis  of  data  from  the  British  station  GAM  near 
Garm,  USSR  for  Soviet  nuclear  explosions 

This  contribution  is  a  follow-up  to  earlier  work  (Ringdal  and  Marshall,  1989; 
Hansen  et  al,  1989;  and  Hansen  and  Ringdal,  1989)  aimed  at  evaluating  the 
stability  of  seismic  Lg  magnitudes  for  yield  estimation  purposes.  In  particu¬ 
lar,  these  efforts  have  involved  analyzing  available  Lg  data  from  Soviet  nuclear 
explosions  at  the  Shagan  River,  Semipalatinsk  test  site,  and  conducting  com¬ 
parative  analyses  of  Lg  and  P  recordings  at  various  seismograph  stations. 

Hansen  et  al  (1989)  analyzed  data  recorded  at  four  digital  stations  installed 
by  IRIS  in  the  Soviet  Union,  and  found  an  excellent  correspondence  between 
Lg  measurements  at  these  stations  and  the  NORSAR  M(Lg)  estimates  pub¬ 
lished  by  Ringdal  and  Marshall  (1989).  Furthermore,  they  noted  the  very  high 
Lg  signal-to-noise  ratio  observed  at  the  IRIS  stations,  in  particular  ARU  and 
GAR,  and  concluded  that  reliable  Lg  measurements  at  these  stations  would  be 
possible  for  explosions  as  small  as  ms  =  4.0,  assuming  normal  noise  conditions. 

Hansen  and  Ringdal  (1989)  extended  the  analysis  to  data  from  the  China 
Digital  Seismograph  Network  (CDSN),  which  is  operated  by  the  USGS  in 
cooperation  with  the  State  Seismological  Bureau,  Beijing.  Two  of  the  CDSN 
stations,  WMQ  in  Urumqi  and  HIA  in  Hailar,  have  particularly  good  Lg  prop¬ 
agation  paths  from  Semipalatinsk,  and  they  based  their  analysis  on  data  from 
these  two  stations. 

In  this  paper,  we  extend  the  analysis  to  data  from  a  broad-band  seismic 
station,  GAM,  installed  very  near  the  IRIS  station  near  Garm,  USSR  (The 
BSVRP  Working  Group,  (1989)).  This  data  supplements  the  previously  sparse 
data  from  GAR  and  allows  the  comparison  of  two  closely  separated  seismic 
stations. 

Fig.  7.2.1  shows  the  locations  of  several  stations  in  the  USSR  and  China 
in  relation  to  the  test  site,  as  weli  as  locations  of  the  NORSAR.  (The  GAM 
station  and  GAR  IRIS  station  are  located  at  the  same  place  on  the  map 
and  indicated  only  by  the  GAR  symbol.)  The  station  GAM  at  a  distance  of 
about  1380  km  shows  excellent  Lg  recordings  of  Semipalatinsk  explosions,  as 
illustrated  by  the  examples  in  Fig.  7.2.2. 

In  the  analysis  of  GAM  Lg  recordings,  we  have  employed  the  exact  same 
procedure  as  described  for  IRIS  data  by  Hansen  et  al  (1989),  and  the  details 
will  not  be  repeated  here.  Data  from  a  total  of  6  Shagan  River  explosions, 
dating  back  to  1988,  were  provided  to  us  for  this  analysis  by  the  BSVRP  Group 
in  Britain.  Table  7.2.1  lists  these  events  along  with  the  estimated  parameters. 

Fig.  7.2.3  shows  a  comparison  of  GAM  and  NORSAR  log  RMS  (Lg)  es¬ 
timates  for  these  6  events.  The  slope  of  the  plot  is  0.92,  and  the  orthogonal 
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standard  deviation  of  the  differences  between  the  two  stations  is  only  0.035 
units.  This  is  essentially  the  same  scatter  found  earlier  by  Hansen  et  al  (1989) 
when  comparing  data  from  NORSAR  and  the  Soviet  station  ARU,  and  con¬ 
firms  the  excellent  stability  of  the  RMS  Lg  estimates. 

As  a  contrast  to  these  well  recorded  events,  Fig.  7.2.4  illustrates  the  ca¬ 
pabilities  of  the  GAM  station  to  record  an  ms(P)  3.8  event  from  the  Shagan 
River  test  site  on  day  270  (September  26)  of  1988.  (This  magnitude  is  based 
on  the  NORSAR  mj(P)  of  4.3  with  an  assumed  regional  correction  of  0.5  units 
for  comparison  to  world  wide  mi,  estimates  and  therefore  must  be  considered 
somewhat  uncertain).  The  unfiltered  broad  band  trace  at  GAM  essentially 
shows  no  signal  for  this  event,  however  the  band  pass  filtered  trace  clearly 
shows  energy  arriving  that  can  be  identified  as  Lg  with  a  signal  to  noise  ratio 
of  about  2.  (Similar  SNR  was  obtained  by  Hansen  ef  al  (1989)  for  the  record¬ 
ing  at  ARU  for  this  event.)  This  SNR  is  near  the  lower  limit  of  about  1.5  for 
allowing  reliable  RMS  Lg  estimates  at  a  single  site. 

Fig.  7.2.5  illustrates  the  stability  of  the  RMS  Lg  amplitudes  by  comparing 
GAM  and  ARU.  These  stations  are  chosen  as  they  are  the  only  pair  for  which 
we  have  Lg  recordings  of  the  m*(P)  3.8  event  shown  in  Fig.  7.2.4  and  so 
illustrate  the  stability  of  measurement  covering  a  span  of  two  full  magnitude 
units.  Here  we  again  have  a  slope  of  very  nearly  one  still  with  an  orthogonal 
standard  deviation  of  only  0.026  logarithmic  units  (i.e.  magnitude  units). 

Fig.  7.2.6  compares  the  signal- to- noise  ratios  (SNR)  (defined  as  RMS  Lg 
signal  to  pre-P  RMS  noise  in  the  0.6  to  3.0  Hz  band)  for  stations  at  various 
distances,  using  5  large  explosions.  The  range  in  magnitude  (ms)  is  from  5.2 
for  the  event  on  day  317  of  1988  to  6.1  for  the  JVE  event  on  day  258  of  1988. 
The  event  on  day  317  indicates  the  minimum  for  which  RMS  Lg  was  measured 
at  NORSAR  at  a  distance  of  about  4200  km  with  a  signal  to  noise  ratio  of 
about  1.1.  For  this  same  event  a  signal  to  noise  ratio  of  about  30  is  observable 
at  ARU  and  GAR  at  a  distance  of  about  1500  km  and  about  80  at  WMQ 
at  a  distance  of  950  km.  Again,  the  event  at  day  258  of  1988  in  Fig.  7.2.6 
(shown  with  the  open  circle  around  a  plus  sign)  shows  an  SNR  gain  of  nearly 
100  between  NORSAR  with  an  SNR  of  3.5  and  WMQ  with  an  SNR  of  331. 
(It  should  be  noted  that  the  low  SNR  for  this  event  at  ARU  is  due  to  the 
fact  that  this  event  was  only  recorded  on  the  low  gain  channel  which  does  not 
adequately  resolve  the  background  noise.)  It  can  be  seen  that  the  SNR  for 
GAM  fits  nicely  to  the  trend  as  a  function  of  distance,  and  actually  is  slightly 
better  than  for  GAR  for  all  common  events. 

In  conclusion,  our  studies  confirm  that  Lg  magnitude  estimates  of  Semi- 
palatinsk  explosions  are  remarkably  consistent  between  stations  widely  dis¬ 
tributed  in  epicentral  distance  and  azimuth.  It  thus  appears  that  a  single  sta¬ 
tion  with  good  signal- to-noise  ratio  can  provide  mt(Lg)  measurements  with  an 
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accuracy  (one  standard  deviation)  of  about  0.03  magnitude  units.  Therefore, 
Lg  signals  appear  to  provide  an  excellent  basis  for  supplying  estimates  of  the 
yields  of  nuclear  explosions  even  down  to  below  one  kiloton,  when  such  sig¬ 
nals  are  recorded  at  high-quality  digital  in-country  seismic  stations,  and  when 
calibrated  by  access  to  independent  (non-seismic)  yield  information  for  a  few 
nuclear  explosions  at  the  test  sites  of  interest.  For  a  review  of  previous  studies 
of  Lg  amplitudes  and  a  more  detailed  account  of  this  work  see  Hansen  et  al 
(1990). 


R.A.  Hansen 
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No. 

Date 

GAM  Lg 

1 

88258 

6.03 

3.184 

2 

88270 

3.8 

1.196 

3 

88317 

5.20 

2.521 

4 

88352 

5.80 

3.034 

5 

89022 

6.0 

3.161 

6 

89043 

5.90 

2.923 

Table  7.2.1  Magnitudes  (mt)  and  log  RMS  Lg  values  at  GAM  for  6  explosions 
analyzed  in  this  study. 


LONGITUDE  (DEG  E) 


Fig.  7.2.1  Map  indicating  the  locations  of  the  Shagan  River  Test  Site,  the 
IRIS  and  British  stations  in  the  USSR,  the  NORSAR  array  in  Norway  and  the 
stations  WMQ  and  HIA  in  China.  The  NORESS  array  is  collocated  near  the 
NORSAR  array,  and  the  station  GAM  is  collocated  near  the  GAR  station. 
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Fig.  7.2.2  Example  of  recordings  from  a  Soviet  nuclear  explosion  (14  Sept 
1988)  at  the  station  GAM.  For  each  of  the  three  components  we  show  the 
unfiltered  trace  (bottom),  the  filtered  trace  (0.6-3.0  Hz)  and  the  120-second 
window  RMS  measure  (top)  as  a  function  of  time. 


Z  Component  RMS  Lg  Comparison 
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Fig.  7.2.3  Comparison  of  log  RMS  Lg  measurements  obtained  at  GAM  and 
NORSAR.  The  standard  deviation  of  the  differences  is  0.035  orthogonal  to  the 
line.  The  dotted  lines  correspond  to  plus  or  minus  two  standard  deviations. 
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Fig.  7.2.4  The  GAM  vertical  component  seismogram  from  the  m /,  3.8  explo¬ 
sion  on  September  26,  1988.  The  lower  trace  is  the  unfiltered  seismogram,  the 
middle  trace  is  the  band  pass  filtered  seismogram  between  0.6  Hz  and  3.0  Hz, 
and  the  upper  trace  is  the  RMS  amplitude  as  a  function  of  time. 


Z  Component  RMS  Lg  Comparison 

S=1.04  1=  -1.10  S.D.=  0.026  N=  6 


Fig.  7.2.5  Comparison  of  log  RMS  Lg  measurements  at  ARU  and  GAM.  The 
slope  of  the  line  is  1.04  and  the  standard  deviation  of  the  misfit  of  the  line  to 
the  data  is  0.026  orthogonal  to  the  line.  The  dotted  lines  correspond  to  plus  or 
minus  two  standard  deviations.  Note  the  remarkable  stability  of  measurement 
between  the  two  stations  over  two  full  magnitude  units. 
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Fig.  7.2.6  Graph  showing  the  variation  of  the  signal-to-noise  ratios  (log  RMS 
minus  log  RMS  noise)  among  GAM,  the  four  IRIS  stations,  the  NORSAR 
array  and  the  CDSN  stations  WMQ  and  HIA.  Epicentral  distance  to  the  test 
Bite  is  plotted  along  the  horizontal  axi*. 
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7.3  Wavefleld  decomposition  for  three-component  seismograms 

Seismograms  at  regional  ranges  are  built  up  from  a  complex  mixture  of  seismic 
propagation  phenomena  involving  multiple  reflections  and  guided  waves  in  the 
crust  (e.g.  the  phases  Pg,  Lg)  and  waves  returned  from  the  uppermost  mantle 
(Pn,  Sn).  In  order  to  detect  and  locate  events  as  well  as  to  assess  the  nature  of 
the  seismic  source  and  propagation  path,  it  is  necessary  to  be  able  to  recognize 
and  characterize  different  parts  of  the  seismic  wavefleld.  Both  phased  array 
techniques  for  single  component  sensors  and  vectorial  analysis  of  3-component 
recordings  can  provide  estimates  of  the  azimuth  and  slowness  of  seismic  phases. 
Although  it  has  been  shown  that  the  accuracy  of  the  3-component  estimates 
of  slowness  and  azimuth  is  generally  lower  than  from  an  array,  especially  for 
S  waves,  a  combination  of  these  approaches  provides  a  more  powerful  tool  to 
estimate  the  propagation  characteristics  of  different  seismic  phases  at  regional 
distances. 

Most  methods  of  using  three-component  data  rely  on  polarization  param¬ 
eters  as  their  main  device  for  characterizing  different  features  on  a  seismic 
wavetrain.  Recently  Jepsen  and  Kennett  (1990)  have  shown  how  it  is  possible 
to  extract  estimates  of  the  relative  amplitudes  of  the  three  incident  wavetypes 
(P,  SV,  SH)  from  three-component  records  under  the  assumption  that  the 
dominant  arrival  at  any  given  time  is  a  single  plane  wave  described  by  its 
slowness  and  azimuth.  This  procedure  depends  on  modelling  the  interaction 
of  the  wavefleld  with  near-receiver  structure.  The  main  contribution,  at  hard 
rock  sites  such  as  NORESS  and  ARCESS,  comes  from  the  interaction  of  the 
wavefleld  with  the  free  surface. 

If  the  slowness  and  azimuth  of  the  incoming  wavefleld  are  known,  then  the 
free  surface  effect  on  an  incident  plane  wave  is  frequency  independent,  but 
may  involve  phase  shifts.  The  extraction  of  the  wavefleld  components  requires 
the  inversion  of  the  matrix  of  interaction  imposed  by  the  free  surface  together 
with  rotation  in  a  horizontal  plane.  For  a  three-component  station  situated  at 
an  array  such  as  ARCESS,  beamforming  over  the  vertical  component  sensors 
can  be  used  to  estimate  the  slowness  and  azimuth  as  a  function  of  time.  These 
array  beam  parameters  may  then  be  used  to  produce  estimates  of  the  P,  SV 
and  SH  contributions  to  the  wavefleld  as  a  function  of  time.  This  process  is 
illustrated  in  Fig.  7.3.1  for  the  3-component  station  C7  at  the  ARCESS  array 
for  an  event  close  to  the  array.  The  array  beam  estimates  were  generated  by 
using  the  broadband  f-k  procedure  on  a  sliding  2  sec  window. 

This  decomposition  of  the  seismic  wavefleld  by  wavetype  as  a  function  of 
time  not  only  has  considerable  benefits  for  the  recognition  of  seismic  phases, 
but  also  provides  a  domain  in  which  the  relative  proportions  of  P,  SV,  and  SH 
can  be  compared  directly,  because  free-surface  amplification  effects  have  been 
removed.  This  information  on  the  current  proportions  of  different  wavetypes 
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summarizes  much  of  the  propagation  processes  between  source  and  receiver 
and  therefore  can  help  to  provide  specific  measures  of  wavefield  character 
which  can  be  beneficial  in  attempts  to  discriminate  between  different  source 
types. 

For  an  individual  3-component  station,  the  estimation  of  apparent  azimuth 
is  usually  more  reliable  than  slowness.  However,  for  regional  phases,  it  proves 
quite  effective  to  use  the  wavefield  decomposition  procedure  with  a  number  of 
fixed  slownesses  designed  to  enhance  different  features  of  the  wavefield  arriving 
at  different  times.  A  slowness  of  0.12  sec/km  gives  very  good  resolution  of  the 
early  P  wave  energy  (see  Fig.  7.3.2)  while  suppressing  later  SV  energy.  P  is 
usually  accompanied  by  a  small  S  component  which  can  bq  reduced  by  optimal 
estimates  of  slowness. 

An  S  wave  slowness  of  0.22  (see  Fig.  7.3.3)  enables  a  clear  identification 
of  the  SV  contribution  to  Lg.  Even  though  the  wavefield  decomposition  pro¬ 
cedure  removes  the  phase  shift  from  free  surface  reflection  imposed  on  SV  at 
large  slownesses,  there  is  commonly  no  close  correspondence  between  the  SV 
and  SH  components.  Rg  waves  are  evident  by  the  presence  of  coupled  P  and 
S  energy. 

An  interesting  by-product  of  a  high  slowness  decomposition  is  an  amplifica¬ 
tion  of  the  P  arrivals  occurring  in  the  notional  SV  trace.  This  occurs  without 
significant  amplification  of  the  background  noise  and  can  lead  to  a  definite 
improvement  in  the  signal- to- noise  ratio.  The  SV  trace  for  high  slowness  can 
therefore  be  investigated  as  a  possible  detector  for  the  P  onset. 

The  wavefield  decomposition  procedure  depends  on  the  specification  of  the 
near  surface  velocities,  but  fortunately  in  general  the  main  results  are  not  very 
sensitive  to  the  values  of  the  velocities.  The  exception  is  for  slownesses  close  to 
the  reciprocal  of  the  surface  P  wave  velocity  (around  0.17  km/s  at  ARCESS) 
where  the  inverse  of  the  free  surface  response  is  a  rapidly  varying  function  of 
slowness.  Also  an  error  of  around  10°  in  azimuth  can  be  tolerated,  the  main 
effect  is  on  the  SH  component. 

Although  slowness  and  azimuth  estimates  from  the  array  are  most  stable 
for  a  particular  frequency  band  (typically  3-5  Hz),  the  f-k  estimates  can  be 
applied  quite  successfully  to  unfiltered  data  or  to  particular  filtered  compo¬ 
nents. 
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Further  examination  of  this  class  of  3-component  analysis  procedures  will 
involve  the  implementation  and  testing  of  the  classification  schemes  of  Jepsen 
and  Kennett  (1990).  Such  schemes  categorize  the  stability,  rectilinearity,  pla¬ 
narity,  etc.  of  the  wavefield  and  use  a  set  of  rules  for  classifying  a  particular 
wavetype  (i.e.  elliptical  P  and  S,  rectilinear  P  and  S,  rectilinear  SH,  Rayleigh, 
etc.). 


B.L.N.  Kennett  &  D.C.  Jepsen,  Research  School  of  Earth  Sciences 
Australian  National  University 
R.A.  Hansen,  NORSAR 
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Fig.  7.3.1.  This  figure  illustrates  the  decomposition  of  a  three  component 
seismogram  into  the  relative  components  of  P,  SV  and  SH.  The  top  three 
traces  are  the  original  three-component  time  series,  the  next  three  traces  are 
for  P,  SV,  and  SH  respectively,  and  the  bottom  two  traces  show  the  values,  as 
a  function  of  time,  of  slowness  and  azimuth  used  in  computing  the  P,  f  and 
SH  contributions.  The  slowness  and  azimuth  were  computed  from  a  2  second 
sliding  window  using  the  entire  array  of  vertical  sensors  and  a  broad  band  f-k 
analysis. 
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Fig.  7.3.2.  This  figure  illustrates  the  decomposition  of  a  three  component 
seismogram  into  the  relative  components  of  P,  SV  and  SH  with  a  fixed  slowness 
of  .12  sec/km  and  azimuth  of  92.  The  top  three  traces  are  the  origin  J  three- 
component  time  series  and  the  bottom  three  traces  are  for  P,  SV,  and  SH 
respectively. 
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Fig.  7.3.3.  This  figure  illustrates  the  decomposition  of  a  three  component 
seismogram  into  the  relative  components  of  P,  SV  and  SH  with  a  fixed  slowness 
of  .22  sec/km  and  azimuth  of  92.  The  top  three  traces  are  the  original  three- 
component  time  series  and  the  bottom  three  traces  are  for  P,  SV,  and  SH 
respectively. 
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7.4  Ray-based  interpretation  of  Lg  azimuth  anomalies  at  NORESS 
Introduction 

In  a  recent  paper,  Bostock  and  Kennett  (1990)  introduced  a  method  for  pre¬ 
dicting  the  propagation  characteristics  of  the  Lg  phase  in  areas  of  complex 
crustal  structure  in  a  semiquantitative  manner.  The  method  relies  on  the 
interpretation  of  Lg  as  constructively  interfering  S  waves  multiply  reflected 
within  the  earth’s  crust.  A  set  of  rays  is  traced  outwards  from  a  point  source 
within  a  crustal  waveguide  of  variable  thickness.  The  geometrical  characteris¬ 
tics  of  the  rays  are  modified  upon  reflection  from  the  crust-mantle  boundary 
and  the  free  surface  according  to  Snell’s  law,  and  can  then  be  monitored  in 
plan  view  by  the  horizontal  projection  of  the  rays  and  the  location  of  mantle 
reflection  points.  This  approach  allows  a  qualitative  description  of  a  variety  of 
complex  propagation  processes  including  mode  coupling,  wavetype  conversion, 
and  lateral  waveguiding  for  situations  which,  because  of  difficulties  in  handling 
boundary  conditions,  cannot  be  adequately  described  by  current  quantitative, 
modal  descriptions  of  the  Lg  phase  (see,  e.g.,  Kennett,  1989).  More  specifi¬ 
cally,  it  is  possible  to  associate  changes  in  the  separation  of  S-wave  reflection 
points  with  a  change  in  the  character  of  the  wavefield  viewed  as  a  sum  of 
higher  mode  surface  waves  since  the  angle  a  ray  makes  with  the  vertical  may 
be  related  to  an  equivalent  phase  velocity.  In  addition,  the  polarization  of 
the  S  waves  may  be  used  as  a  measure  of  the  conversion  between  Love  and 
Rayleigh  waves. 

Although  the  method  has  been  used  to  explain  the  observed  patterns  of 
Lg  propagation  in  central  Asia  (Bostock  and  Kennett,  1990)  and  the  south¬ 
western  United  States  (Kennett  et  al,  1990),  its  success  is  highly  dependent 
upon  the  knowledge  of  crustal  model.  In  these  previous  studies,  variation  in 
crustal  thickness  was  considered  the  dominant  factor  influencing  the  character 
of  Lg  on  seismograms,  and  crustal  models  based  on  varying  degrees  of  iso¬ 
static  compensation  were  adopted  in  accordance  with  published  geophysical 
and  tectonic  evidence.  However,  there  are  a  number  of  other  factors  which 
may  become  significant  in  general  applications.  For  instance,  in  areas  of  thick 
sedimentary  cover  it  may  be  advisable  to  consider  that  portion  of  the  crust 
comprising  basement  alone.  In  general,  the  effects  of  vertical  velocity  struc¬ 
ture  and  systematic  lateral  variations  therein  will  generally  be  insignificant 
at  all  but  the  lowest  phase  velocities,  however  abrupt  lateral  heterogeneity 
in  physical  properties  will  undoubtedly  play  an  important  role  in  distorting 
the  observed  wavefield.  A  truly  quantitative  representation  of  Lg  propagation 
using  the  ray  method  is  only  possible  when  comprehensive  information  on  all 
contributing  factors  is  incorporated  in  the  construction  of  our  crustal  models. 
Notwithstanding,  in  many  cases  valuable  insight  may  be  gained  by  considering 
simply  the  variation  in  crustal  thickness,  in  which  case  an  accurate  knowledge 
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of  the  depth  to  Moho  is  essential  since  fluctuations  in  this  lower  boundary  are, 
in  general,  more  pronounced  than  those  at  the  free  surface. 

Scandinavia  crustal  model 

Our  objective  in  this  study  is  to  employ  the  ray  method  as  a  tool  in  assessing 
the  effects  of  gross  crustal  structure  on  Lg  propagation  in  Scandinavia  and 
adjacent  areas  to  the  east,  and  relating  the  predicted  behavior  to  azimuth 
anomalies  observed  at  the  small-aperture  NORESS  array.  We  will  therefore 
consider  a  single-layer  crustal  model  over  the  area  0°-35°E,  550-75°N  com¬ 
prising  all  of  Scandinavia  and  Finland  and  parts  of  the  Soviet  Baltic  states. 
The  modelled  free  surface  boundary  is  a  smoothed  (81  point  average)  version 
of  relief  data  from  the  ET0P05  world  topography  data  base  at  5’  intervals. 
The  crust-mantle  boundary  was  constructed  by  interpolating  the  Moho  map 
of  Ruud  (pers.  comm.)  which  combines  information  from  a  variety  of  geo¬ 
physical  studies  (seismic  reflection,  refraction,  gravity,  etc.)  and  is  thought  to 
be  an  accurate  representation  of  the  main  low-wavelength  structure  over  the 
region.  The  main  topographic  feature  is  the  Caledonide  Range  paralleling  the 
coast  of  Norway,  which  rises  to  an  elevation  of  2000+  m  over  a  considerable 
area.  The  Moho  map  shown  in  Fig.  7.4.1  presents  a  slightly  more  compli¬ 
cated  picture.  The  mean  sea  level  depth  to  Moho  decreases  rapidly  off  the 
Norwegian  coast  to  typical  values  for  oceanic  crust  (~  10  km).  Interestingly, 
the  Caledonides  are  not  mirrored  by  any  significant  root  on  the  Moho;  rather, 
crustal  thicknesses  generally  increase  steadily  as  one  approaches  the  center  of 
the  Baltic  Shield.  There  are  several  areas  of  increased  crustal  thickness  (> 
50  km);  notably  along  the  central  and  northern  portions  of  the  Soviet-Finnish 
border,  and  a  saddle-shaped  structure  straddling  the  eastern  coast  of  Sweden. 
Moho  depths  tend  to  decrease  further  south  and  east  approaching  the  Russian 
platform. 

Ray  analysis 

Mykkeltveit  et  al  (1989)  have  assembled  a  data  base  of  Lg  azimuth  anomalies 
observed  at  NORESS  from  a  variety  of  sources  at  distances  from  100  to  1000 
km  which  indicates  that  in  some  regions  there  is  extreme  variability  in  the 
magnitude  and  sign  of  azimuth  anomalies  with  even  minor  changes  in  source 
position.  We  will  examine  possible  causes  for  this  anomaly  distribution  using 
the  ray  technique  discussed  above. 

Perhaps  the  most  interesting  region  in  terms  of  anomalous  behavior  is 
located  along  the  southern  coast  of  Finland  on  the  Gulf  of  Bothnia.  Here 
Lg  waves  from  sources  separated  by  distances  as  small  as  2°-3c  in  latitude 
exhibit  strikingly  different  azimuthal  anomalies.  To  investigate  this  behavior, 
we  present  ray  diagrams  generated  for  i)  different  phase  velocities  at  a  single 
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source  location  (21°E,  60°N)  in  Fig.  7.4.2,  and  ii)  different  source  locations  in 
this  general  area  (21°E,  59.0°-62.0“N)  at  constant  phase  velocity  (4.0  km/s) 
in  Fig.  7.4.3. 

We  first  describe  the  variable  phase  velocity  diagrams  shown  in  Fig.  7.4.2. 
For  lower  phase  velocities  (e.g.,  3.8  km/s  —  Fig.  7.4.2a)  guided  wave  prop¬ 
agation  is  generally  less  affected  by  lateral  variation  in  crustal  thickness  as 
an  individual  ray  undergoes  fewer  bounces  over  a  given  distance.  The  fig¬ 
ure  also  indicates  that  a  large  proportion  of  the  rays  undergo  reflection  at 
the  continental  margin  off  western  Norway  accompanied  by  conversion  to  Sn 
mantle  phases  into  the  oceanic  crust  (denoted  by  solid  diamonds).  Of  specific 
relevance  is  the  saddle-shaped  low  in  Moho  relief  centered  roughly  midway 
between  NORESS  and  the  source,  and  occupying  a  significant  portion  of  the 
total  path  (see  also  Fig.  7.4.1).  Note  that  the  axis  of  the  ridge  separating  the 
two  pockets  is  roughly  colinear  with  a  line  joining  NORESS  and  the  source. 
We  recall  that  zones  of  increased  crustal  thickness  tend  to  behave  as  attrac¬ 
tors,  pulling  rays  inward  in  much  the  same  fashion  as  high  velocity  zones 
in  conventional  body  wave  raytracing.  Hence  as  a  general  observation,  rays 
launched  at  all  3  phase  velocities  tend  to  be  drawn  away  from  the  Moho  ridge 
leaving  windows  through  which  fewer  rays  pass.  These  windows  overlap  for 
the  three  ray  diagrams  (Figs.  7.4.2a,  b,  c)  but  do  not  exactly  coincide  since 
the  phase  velocity  dictates  the  exact  location  of  the  basement  reflection  points 
for  a  given  ray  and  their  resulting  ray  curvature.  The  position  of  NORESS 
relative  to  these  alternate  low  and  high  ray  density  windows  is  therefore  a 
function  of  phase  velocity,  and  we  might  accordingly  expect  the  azimuth  of 
the  observed  wavetrain  to  vary  considerably  with  time  as  the  individual  phase 
velocity  components  arrive  at  different  group  velocities.  In  addition,  observed 
azimuths  are  likely  to  depend  significantly  on  the  frequency  characteristics  of 
the  particular  source-excitation  function. 

'  ve  now  consider  the  effect  of  shifting  the  source  by  1°  increments  along 
the  21°  E  meridian  on  ray  diagrams  generated  at  a  fixed  phase  velocity  of  4.0 
km/s  as  shown  in  Fig.  7.4.3.  Note  that  even  small  changes  in  source  position 
can  significantly  alter  the  source-heterogeneity-receiver  geometry.  The  ray 
diagrams  reflect  this  fact  by  exhibiting  considerable  variation  in  ray  density 
west  of  the  main  heterogeneity:  the  saddle-shaped  low  in  Moho  relief  behaves 
as  a  lens  which  focuses  and  defocuses  ray  bundles  in  a  way  that  is  strongly 
dependent  on  the  location  of  the  illuminating  source. 

Finally  we  examine  a  more  distant  source  at  the  head  of  the  Gulf  of  Estonia. 
The  ray  diagram  in  Fig.  7.4.4  indicates  that  the  absence  of  major  heterogeniety 
along  the  initial  part  of  the  ray  paths  results  in  a  ray  pattern  which  remains 
quite  coherent  into  eastern  Sweden.  Rays  in  the  vicinity  of  the  heterogeneity 
are  more  nearly  parallel  than  those  from  closer  sources  (as  in  Figs.  7.4.2  and 
7.4.3)  and  tend  not  to  diverge  as  markedly  upon  transmission.  Nevertheless, 


the  effects  of  the  Moho  saddle  are  still  manifest  in  a  focusing  and  defocusing 
of  rays  which,  for  the  particular  phase  velocity  employed,  leaves  NORESS  in 
a  region  of  low  ray  density. 

Discussion 

Initial  attempts  to  relate  ray-calculated  and  observed  azimuth  anomalies  yield¬ 
ed  somewhat  inconsistent  results,  but  upon  careful  examination  of  the  ray 
behavior  it  becomes  evident  that  the  nature  of  the  Moho  heterogeneity  in 
the  vicinity  of  NORESS  is  probably  largely  responsible.  In  areas  of  broader 
scale  heterogeneity,  ray  diagrams  are  characterized  by  more  systematic  and 
gradual  changes  in  ray  pattern  with  respect  to  both  source  location  and  phase 
velocity  (see  for  example  Bostock  and  Kennett,  1990).  The  complexity  of  the 
Moho  saddle  east  of  NORESS,  its  comparatively  complicated  geometry  over  a 
region  of  restricted  spatial  extent,  results  however  in  ray  diagrams  which  are 
considerably  more  sensitive  to  these  two  parameters.  The  observations  made 
above  would  appear  to  bear  important  implications  to  the  analysis  of  crust 
and  upper  mantle  phases  originating  from  parts  east  of  the  array,  especially  for 
more  proximal  sources.  The  study  suggests  that  the  Moho  saddle,  a  feature 
which  appears  to  be  well  defined  and  documented  in  a  number  of  independent 
investigations,  will  play  a  significant  role  in  distorting  the  wavefield  and  is 
probably  responsible  in  large  part  for  the  Lg  azimuth  anomalies  observed  at 
NORESS  from  sources  to  the  east. 


M.G.  Bostock  fe  B.L.N.  Kennett,  Research  School  of  Earth  Sciences 
Australian  National  University,  Canberra 
S.  Mykkeltveit,  NORSAR 
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Fig.  7.4.1.  Moho  relief  over  Scandinavia,  Finland  and  the  Soviet  Baltic  states 
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Fig.  7.4.2.  Ray  diagrams  for  a  source  located  at  210E,  60°N  at  three  phase 
velocities  a)  3.8  km/s,  b)  4.0  km/s  and  c)  4.2  km/s.  Note  alternating  high 
and  low  ray-density  windows  resulting  from  Moho  saddle  east  of  NORESS. 
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Fig.  7.4.3.  Ray  diagrams  for  3  sources  located  at  a)  21°E,  59°N,  b)  21°E, 
61°N  and  c)  21°E,  62°N  at  a  phase  velocity  of  4.0  km/s.  Note  variability  of 
ray  patterns  with  source  location. 
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Gulf  of  Estonia  2 
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Fig.  7.4.4.  Ray  diagram  for  a  source  in  the  eastern  Gulf  of  Estonia  at  a  phase 
velocity  of  4.0  km/s.  Note  the  relative  coherence  in  the  ray  patte'.n  east  of 
the  Swedish  coast. 
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7.5  Earthquake  reporting  capabilities  in  Fennoscandia  as  in¬ 
ferred  from  IAS  data 

For  regional  seismic  arrays  and  networks  in  populated  areas,  operating  at 
frequencies  above  1  Hz,  it  is  commonly  observed  that  detections  become  in¬ 
creasingly  dominated  by  man-made  events  as  the  signal- to-noise  ratio  (SNR) 
is  decreased.  It  is  also  known  that  the  spatio-temporal  characteristics  of  nat¬ 
ural  earthquakes  are  often  quite  different  from  man-made  events.  Most  of  the 
events  of  antropogenic  origin  are  chemical  explosions  of  some  sort. 

The  Intelligent  Array  System  (IAS)  installed  at  NORSAR  is  in  its  present 
version  a  system  for  routine  processing  and  analysis  of  data  from  the  regional 
arrays  NORESS  and  ARCESS  in  Norway  (Bache  et  al,  1990;  Bratt  et  al ,  1990). 
From  September  29,  1989,  and  to  the  end  of  the  year  the  system  was  operated 
at  the  Center  for  Seismic  Studies  (CSS)  in  Arlington,  Virginia,  USA,  and  from 
the  beginning  of  1990  at  NORSAR. 

The  purpose  of  the  present  study  is  to  evaluate  the  data  from  the  first 
six  months  of  operation  of  the  IAS,  from  the  point  of  view  of  earthquake 
identification.  In  doing  this,  we  use  available  information  about  the  seismicity 
of  the  region  (Bungum  et  al,  1990),  but  with  no  prior  assumptions  otherwise 
about  where  to  expect  the  man-made  event  . 

Data  and  data  analysis 

The  IAS  system  is  fully  automatic,  but  also  supplemented  by  a  thorough 
analyst  review  and  evaluation.  During  the  6  months  analyzed  here,  a  total  of 
3813  located  events  were  accepted  by  the  analysts,  who  refined  in  one  way  or 
another  the  solutions  for  all  but  669  of  these.  To  further  improve  the  quality 
of  the  event  locations  selected  for  the  present  analysis,  we  have  accepted  only 
solutions  where  the  standard  error  of  one  observation  (sdobs)  is  less  than  3 
seconds  and  where  the  standard  error  in  time  (stt)  is  less  than  10  seconds 
(Bratt  and  Bache,  1988).  This  reduced  the  number  of  events  to  2454,  of  which 
all  but  9  are  located  within  the  greater  Fennoscandian  region  of  54-76°N, 
0-40°E. 

The  spatial  distribution  of  these  events  is  shown  in  Figs.  7.5.1a-c,  revealing 
a  pattern  quite  different  from  the  known  seismicity  in  the  area,  and  with  a 
non-uniform  tirae-of-day  (TOD)  distribution  as  shown  in  Fig.  7.5. Id.  This 
distribution  contrasts  strongly  to  that  expected  for  earthquakes,  with  no  de¬ 
pendence  on  time-of-day.  There  are  two  prominent  peaks  in  the  distribution, 
right  ahter  noon  and  around  midnight.  In  fact,  the  amount  of  man-made 
events  (explosions)  in  the  data  base  is  larger  than  what  is  indicated  by  Fig. 
7.5.1d,  since  different  areas  should  be  expected  to  have  their  explosions  peaks 
at  different  times  of  the  day,  thereby  cancelling  each  other  to  some  extent. 
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Our  aim  in  the  present  study  has  been  to  identify  areas  with  a  signifi¬ 
cant  number  of  explosions,  and  to  flag  these  events  in  the  data  base.  Our 
approach  in  this  respect  has  been  first  to  identify  epicentral  areas  where  the 
events  are  concentrated  or  clustered.  A  non-uniform  TOD  distribution  and/or 
a  non-exponential  magnitude- frequency  (MF)  distribution  then  indicates  the 
presence  of  man-made  events  within  that  area.  The  best  way  to  do  this  anal¬ 
ysis  is  to  start  with  large  areas  or  regions,  and  then  to  zoom  in  on  smaller 
areas  later. 

For  areas  thus  identified  to  contain  a  certain  amount  of  explosions,  we 
define  one  or  more  hourly  intervals  within  which  all  events  below  a  certain 
magnitude  should  be  flagged  as  explosions  in  the  data  base.  These  hourly 
intervals  coincide  with  peaks  in  the  TOD  distributions,  while  the  magnitude 
thresholds  depend  on  the  character  of  the  MF  distributions  within  as  well  as 
outside  of  these  intervals.  The  parameters  should  be  tuned  in  such  a  way 
as  to  yield  (after  flagging  of  explosions)  both  TOD  and  MF  distributions  for 
the  remaining  events  that  are  compatible  with  what  should  be  expected  for 
earthquakes. 

After  screening  or  filtering  the  entire  data  base  in  this  way,  we  evaluate 
what  this  filtering  represents  in  terms  of  magnitude  thresholds  for  accepting 
events  (as  natural  earthquakes)  from  particular  areas,  at  different  times  of 
the  day.  These  thresholds  reflect  in  the  present  context  the  capabilities  for 
reporting  natural  earthquakes,  limited  by  a  noise  level  which  in  this  case  is 
determined  by  the  occurrence  of  man-made  events. 

The  actual  analysis  for  the  IAS  data  base  involved  the  testing  of  a  large 
number  of  spatial  windows,  in  order  to  identify  areas  with  non-uniform  TOD 
distributions  and  non-exponential  MF  distributions.  The  result  of  this  analysis 
for  the  present  data  base  of  2454  Fennoscandian  event  locations  are  the  52 
TOD  filters  defined  in  Table  7.5.1.  The  spatial  windows  are  shown  also  in 
Figs.  7.5.1a-c,  and  the  TOD  distributions  for  a  selection  of  these  windows  are 
shown  in  Table  7.5.2. 

In  some  cases,  a  spatial  window  is  defined  entirely  within  the  area  covered 
by  another  window.  In  that  case,  the  smaller  window  has  a  higher  magni¬ 
tude  threshold.  In  other  cases,  one  window  has  two  different  hourly  intervals 
connected  to  it,  which  is  done  by  defining  two  different  windows.  After  the 
first  47  windows  were  defined  (only  those  are  shown  in  Figs.  7.5.1a-c),  it  was 
found  necessary  to  define  also  four  larger  regional  windows  (No  48-51),  but 
with  magnitude  limits  below  the  smaller  windows  covered  inside.  In  addition, 
the  last  window  (No  52)  covers  the  entire  area  under  analysis,  for  the  purpose 
of  defining  a  lower  magnitude  limit  for  all  the  events  (here  set  to  0.5).  From 
Table  7.5.1  it  is  seen  that  the  other  windows  have  magnitude  limits  ranging 
from  1.5  to  3.3. 
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After  TOD  filtering,  there  are  75  events  remaining  of  the  original  2454, 
or  about  3%.  These  events  are  plotted  in  Fig.  7.5.2,  with  TOD  distributions 
as  shown  in  Fig.  7.5.1d  and  in  in  the  last  column  of  Table  7.5.2.  It  is  likely, 
when  studying  the  TOD  distribution  and  also  when  comparing  with  earlier 
seismicity  maps,  that  some  explosions  still  remain  in  the  data  base.  This 
is  because  they  occur  in  areas  with  few  other  events,  because  they  occur  at 
TODs  shared  by  few  other  events  from  the  same  area,  or  because  they  have 
magnitudes  above  other  events  from  the  same  area  and  with  the  same  TODs. 

While  it  is  likely  that  the  TOD  filtering  has  not  removed  all  of  the  ex¬ 
plosions,  it  is  even  more  certain  that  some  natural  earthquakes  have  been 
removed.  The  latter  is  evident  from  the  fact  that  what  we  have  done  is  simply 
to  remove  all  of  the  events  below  a  certain  threshold  magnitude,  a  threshold 
which  has  been  made  dependent  on  region  and  TOD.  In  the  cases  when  natu¬ 
ral  earthquakes  occur  within  these  spatio-temporal-magnitude  windows,  they 
will  of  course  be  flagged  as  explosions. 

Earthquake  reporting  capabilities 

The  results  of  the  above  analysis  for  six  months  of  available  IAS  data 
determines  a  threshold  for  the  reporting  of  natural  earthquakes,  provided  that 
no  discrimination  analysis  of  individual  events  has  been  performed.  The  key 
to  this  threshold,  and  the  way  in  which  it  varies  in  time  and  space,  is  found 
in  the  definition  of  the  TOD  filters  (Table  7.5.1). 

We  have  developed  a  procedure  for  time  dependent  spatial  contouring  of 
this  threshold  by  defining  first  an  hourly  time  interval  and  a  regular  spatial 
grid.  For  each  grid  point,  we  loop  through  all  TOD  filters  in  order  to  identify 
those  that  contain  the  point,  and  a  hit  is  declaired  if  there  is  some  overlap 
between  the  selected  hourly  interval  and  the  filter’s  hourly  interval.  The  mag¬ 
nitude  threshold  for  the  selected  grid  point  is  then  taken  from  the  filter  with 
the  highest  threshold.  This  procedure  is  repeated  for  all  points  in  the  grid, 
all  of  which  are  associated  with  a  threshold  magnitude  applicable  within  the 
selected  TOD  interval.  A  threshold  contouring  can  then  be  performed  on  this 
basis. 

In  the  present  case,  we  have  defined  a  grid  with  steps  0.5°  in  latitude 
and  1°  in  longitude.  The  resulting  contours  (with  some  spatial  smoothing)  is 
shown  in  Fig.  7.5.2a-b  for  a  time  period  (1030-1430  UTC)  during  the  peak  of 
the  working  hours,  and  in  Fig.  7.5.2c-d  for  a  time  period  (0030-1430  UTC) 
during  the  most  quiet  part  of  the  night.  The  difference  between  the  two  is  very 
significant:  while  the  noise  level  (reporting  threshold)  at  night  time  extends 
up  to  1.5  in  magnitude  only  in  very  few  areas,  the  level  in  the  middle  of  the 
day  extends  up  to  2.5  over  fairly  large  areas. 
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Discussion  and  conclusions 

The  hourly  time  period  used  in  the  examples  above  is  somewhat  arbitrarily 
set  to  4  hours.  The  length  of  the  time  interval  should  reflect  the  particular 
purpose  of  the  analysis,  but  it  should  also  reflect  the  time  resolution  assumed 
in  the  definition  of  the  TOD  filters.  A  similar  resolution  consideration  applies 
also  to  the  determination  of  the  grid  size,  which  should  reflect  the  purpose  of 
the  analysis  as  well  as  the  resolution  assumed  in  the  definition  of  the  TOD 
filters. 

One  of  the  limitations  with  the  present  data  base  is  that  it  is  based  on 
observed  data  from  only  two  regional  arrays.  However,  even  if  this  leads  to  a 
certain  variation  in  both  detection  thresholds  and  location  precisions  over  the 
area  covered  in  this  study,  it  should  not  affect  the  results  significantly.  The 
reason  for  this  is  partly  that  the  poorest  locations  have  been  removed  prior  to 
the  analysis,  but  first  of  all  that  a  substantial  amount  of  the  explosions  are  in 
fact  located  more  or  less  randomly  outside  of  the  main  known  explosion  sites. 
Better  locations  should  therefore  only  be  expected  to  lead  to  minor  changes 
in  the  space-time  organization  of  the  data. 

The  spatial  windows  as  defined  in  this  study  are  very  simple,  reflecting 
to  some  extent  the  limitations  in  location  precisions.  A  natural  refinement 
here  would  be  to  define  windows  in  terms  of  polygones,  an  option  which  in 
fact  is  available  in  the  sorting  program  used  in  this  study.  In  the  present  first 
order  approach,  however,  simple  rectangular  spatial  windows  are  considered 
satisfactory. 

Another  natural  refinement  of  the  window  definitions  would  be  to  include 
weekly  and  seasonal  variations,  both  of  which  are  known  to  be  considerable. 
The  weekly  variation  shows  for  most  of  the  areas  a  stable  weekend  minimum 
(Bratt  et  al,  1990),  while  the  seasonal  variations  probably  are  more  compli¬ 
cated.  An  inclusion  in  the  filter  definitions  of  such  variations  is  quite  straight¬ 
forward,  however. 

The  explosion  filters  and  associated  noise  levels  defined  in  this  work  are  of 
course  applicable  in  general  only  for  the  time  period  for  which  they  has  been 
derived.  However,  only  minor  adjustments  should  normally  be  expected  to 
be  necessary  in  order  to  apply  the  present  results  to  other  time  periods.  It 
would  be  desirable  to  combine  such  an  extention  with  the  refinements  discussed 
above.  Similarly,  the  approach  taken  here  should  also  be  easy  to  apply  to  any 
new  area  from  which  event  location  data  of  similar  nature  are  available. 

What  has  been  shown  here  is  that  the  detectability  now  possible  through 
a  network  of  regional  arrays  cannot  be  fully  utilized  unless  a  substantial  effort 
is  being  directed  into  the  problem  of  event  identification  and  discrimination. 
Being  interested  in  earthquakes  or  explosions  are  two  sides  of  the  same  ques- 
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tion  in  this  respect,  pursuable  through  the  same  approaches,  methods  and 
algorithms. 

To  be  most  useful,  the  discrimination  algorithms,  whether  they  are  based 
on  focal  depth,  spectral  features,  spectral  ratios,  etc.  (e.g.  Hedlin  et  at ,  1990), 
should  to  the  largest  possible  extent  be  included  in  the  routine  analysis  of  the 
events,  with  options  for  refinements  and  improvements  as  part  of  the  analyst 
reviews  as  well  as  in  subsequent  offline  analyses.  These  results  should  follow 
the  events  into  the  data  base  and  be  expressed  in  ways  that  could  facilitate 
the  computation  of  probabilities  for  the  events  being  natural  earthquakes  or 
not. 

In  conclusion,  we  have  found  through  a  regionalized  time-of-day  analysis  of 
six  months  of  IAS  data  from  the  regional  arrays  NORESS  and  ARCESS  shows 
that  about  97%  of  the  more  well-located  events  probably  are  of  man-made 
origin.  Based  on  the  derived  information  about  where  and  when  the  man¬ 
made  events  occur,  contour  maps  (in  magnitude)  of  associated  capabilities 
for  reporting  natural  earthquakes  are  provided,  for  different  times  or  hourly 
intervals  of  the  day.  The  magnitude  limits  vary  from  about  3.0  as  a  maximum 
in  some  mining  areas  and  down  to  less  than  1.0  for  western  and  northern 
Norway  offshore  areas,  where  most  of  the  man-made  disturbances  are  still 
below  the  detection  level. 

These  results  call  for  a  dedicated  effort  into  event  discrimination  work, 
with  algorithms  to  be  included  preferably  already  in  the  routine  analysis  of 
the  data. 
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THble  7.5.1.  Time-of-day  (TOD)  filters  as  developed  in  this  study  for  the 
purpose  of  removing  presumed  explosions  from  the  IAS  data  base  of  event 
locations  in  Fennoscandia.  The  columns  are  filter  (spatio-temporal  window) 
number,  latitude  and  longitude  limits,  magnitude  limit  and  time-of-day  limits. 
The  last  of  the  filters  (52)  covers  the  entire  region  under  analysis,  thereby 
defining  a  lower  magnitude  threshold  for  accepting  any  event. 
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Table  7.5.2.  Time-of-day  (TOD)  distributions  for  the  spatio-temporal  win¬ 
dows  (filters)  defined  in  Table  7.5.1  that  have  more  than  10  events  in  any 
particular  hourly  interval.  The  two  last  columns  (shown  also  if  Fig.  7.5. Id) 
give  the  distribution  in  the  original  data  base  and  in  the  filtered  one,  respec¬ 
tively. 
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Fig.  7.5.2.  Contour  maps  for  the  magnitude  thresholds  defined  by  the  TOD 
filters  developed  in  this  study.  The  four  maps  correspond  to  (a)  magnitude 
level  of  1.8  for  the  hours  between  1030  and  1430  UTC,  (b)  magnitude  2.4  for 
1030-1430,  (c)  magnitude  0.6  for  0030-0430  and  (d)  magnitude  1.7  for  0030- 
0430.  Under  certain  conditions,  these  contours  indicate  the  day-time  reporting 
capabilities  for  natural  earthquakes. 


7.6  Application  of  the  threshold  monitoring  method 

The  concept  of  threshold  monitoring  was  introduced  by  Ringdal  and  Kvaerna 
(1989)  as  a  method  of  monitoring  the  seismic  amplitude  levels  for  the  pur¬ 
pose  of  using  this  information  to  assess  the  largest  size  of  events  that  might 
go  undetected.  In  an  effort  to  demonstrate  the  capabilities  of  this  threshold 
monitoring  concept,  a  preliminary  version  has  been  implemented  into  the  In¬ 
telligent  Monitoring  System  (IMS)  (Bache  et  at,  1990).  A  demonstration  of 
this  implementation  was  given  at  the  Symposium  on  Regional  Seismic  Arrays 
and  Nuclear  Test  Ban  Verification,  held  in  Oslo  in  February  1990.  In  the 
following,  we  will  present  figures  from  that  demonstration,  as  well  as  a  brief 
description  of  the  method. 

Method  description 

The  basic  idea  behind  the  threshold  monitor  is,  for  any  given  point  in  time, 
to  infer  the  upper  magnitude  limit  of  a  possible  seismic  event  at  a  given  geo¬ 
graphical  location.  By  combining  observations  of  the  amplitude  of  the  seismic 
data  at  different  arrays  and/or  single  stations,  we  can  apply  the  formalism 
developed  by  Ringdal  and  Kvaerna  (1989)  to  compute  an  upper  magnitude 
limit  based  on  the  network. 

In  order  to  apply  this  method  the  following  procedure  is  required: 

•  For  each  location-station-phase  combination,  estimate  continuously  the 
seismic  amplitude  levels.  If  the  station  is  an  array,  we  use  STA  values  of 
filtered  beams  to  represent  the  amplitude  levels.  The  steering  parameters 
of  the  beams  will  then  correspond  to  the  apparent  velocity  and  azimuth 
of  the  actual  phase.  The  filter  bands  are  chosen  such  that  good  SNR  is 
ensured.  If  the  observation  unit  is  a  single  station,  the  STA  values  are 
computed  from  a  filtered  channel. 

•  When  considering  a  potential  event  at  a  given  time  and  location,  mea¬ 
sure  the  seismic  amplitude  levels  at  the  expected  arrival  times  for  the 
relevant  seismic  phases.  The  travel  times  for  each  phase  can  be  taken 
from  standard  travel  time  tables,  or  by  processing  events  with  known 
location  and  origin  time. 

•  In  order  to  relate  the  STA  observations  to  actual  magnitude  estimates, 
apply  the  formula 

m  =  log(STA)  +  b 

where  b  is  a  correction  factor  for  each  location-station-phase  combina¬ 
tion.  The  correction  factors  can  be  obtained  by  processing  events  with 
known  magnitudes,  or  by  using  standard  attenuation  values. 
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•  For  assessing  the  significance  of  our  magnitude  estimates,  assume  that 
they  are  sampled  from  a  normal  distribution  with  a  given  standard  de¬ 
viation.  Based  on  experience  with  signal  amplitude  variation  across  the 
NORSAR  array,  we  have  used  a  preliminary  value  of  0.2  for  a  small 
epicentral  area. 

•  The  magnitude  limits  computed  by  this  algorithm  are  tied  to  a  given 
confidence  level,  here  set  to  0.9.  This  means  that  the  estimated  limits 
represent  the  largest  magnitude  of  a  possible  hidden  event,  in  the  sense 
that  there  is  at  least  a  90  per  cent  probability  that  one  or  more  of  the 
observed  amplitude  values  would  be  exceeded  by  the  signals  from  an 
event  with  magnitude  above  these  limits. 

Interfering  events 

To  illustrate  the  capability  of  the  threshold  monitor  in  the  presence  of  an 
interfering  event,  we  assume  a  situation  as  shown  in  Fig.  7.6.1.  A  teleseismic 
event  is  interfering  with  a  local  mining  explosion,  causing  large  amplitudes 
at  the  expected  arrival  time  of  Pn.  We  also  assume  that  we  are  monitoring 
the  source  region  of  the  actual  mining  explosion,  using  appropriate  calibration 
values  for  the  seismic  phases  considered. 

If  we  now  compute  an  upper  magnitude  limit  only  from  the  amplitude  level 
at  the  time  of  the  Pn  arrival  (where  the  interfering  event  is  added),  we  will 
necessarily  conclude  that  a  relatively  large  explosion  may  have  occurred.  If 
we  in  addition  bound  the  event  magnitude  by  the  undisturbed  Lg  amplitude, 
we  will  get  a  much  lower  value. 

An  actual  interfering  event  situation  may  in  fact  be  more  complex  than 
this,  e.g.  one  station  may  be  contaminated  with  high  amplitudes  for  a  long 
time  period,  or  may  even  not  record  data.  In  such  cases,  amplitudes  recorded 
at  other  stations  may  put  strong  constraints  on  the  upper  magnitude  limit, 
and  thereby  exclude  the  possible  occurence  of  a  strong  event. 

Beamforming 

The  beamforming  capability  of  regional  arrays  is  efficiently  exploited  by  the 
threshold  monitor,  both  for  P  and  S-phases.  The  STA  traces  used  in  the 
computation  of  the  upper  magnitude  limits  are  derived  from  the  amplitudes 
of  filtered  steered  beams.  At  NORESS,  we  know  that  forming  Pn  beams  will 
reduce  the  amplitude  level  of  pure  noise  by  about  14  dB  (Kvserna,  1989). 
For  optimally  steered  Pn  beams,  the  signal  loss  will  be  less  than  3  dB  for 
frequencies  below  6  Hz.  This  will  cause  the  calibration  factors  (b-values)  to 
increase  with  the  same  amount,  but  we  will  still  lower  the  upper  magnitude 
limit  by  more  than  10  dB  (0.5  magnitude  units). 
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For  Lg,  we  expect  the  noise  suppression  to  be  somewhat  less  than  for  Pn, 
and  we  also  expect  a  higher  signal  loss.  Nevertheless,  a  reduction  of  the  upper 
magnitude  limit  by  more  than  6  dB  seems  to  be  feasible  for  Lg. 

Signals  arriving  from  events  outside  the  monitoring  region  may  have  their 
amplitudes  significantly  reduced  by  beamforming.  The  amount  of  amplitude 
reduction  depend  on  the  difference  between  the  slowness  vectors  of  the  beam 
and  the  arriving  signal,  as  well  as  on  the  beampattern  of  the  array. 

Implementation 

In  the  present  demonstration  we  use  data  from  the  three  regional  arrays  in 
Fennoscandia:  NORESS,  ARCESS  and  FINESA.  As  target  regions  for  thresh¬ 
old  monitoring,  we  have  chosen  10  active  mining  regions  in  Fennoscandia  and 
Western  USSR,  as  well  as  the  Soviet  nuclear  test  site  at  Shagan  River  (see 
Table  7.6.1). 

For  the  10  mining  regions,  the  phases  Pn  and  Lg  at  each  array  are  being 
used  to  infer  the  upper  magnitude  limits.  Calibration  factors  are  found  by 
processing  events  with  magnitudes  reported  in  the  regional  seismic  bulletin 
of  the  University  of  Helsinki.  Travel  times  and  beamsteering  parameters  are 
obtained  from  the  same  processing.  The  filter  bands  for  Pn  and  Lg  are  re¬ 
spectively  3-5  Hz  and  1.5-3.5  Hz,  and  the  STA  sampling  interval  is  set  to  1 
second.  For  Pn,  the  STA  integration  window  is  2  seconds,  and  maximum  STA 
is  chosen  within  ±2  seconds  of  the  predicted  arrival  time.  The  reason  why 
the  amplitude  level  is  represented  by  the  maximum  STA  within  a  certain  time 
tolerance,  is  that  each  target  point  represents  a  finite  region,  e.g.  10x10  km, 
and  that  phases  from  events  occuring  outside  the  center  point  will  have  some¬ 
what  different  travel  times.  The  STA  integration  window  for  Lg  is  10  seconds, 
with  maximum  chosen  within  ±3  seconds  of  the  predicted  arrival  time.  For 
distances  less  than  500  km,  however,  the  integration  window  for  Lg  is  reduced 
to  5  seconds. 

For  the  Shagan  River  test  site,  calibration  factors,  travel  times  and  beam¬ 
steering  parameters  for  the  phases  P,  PP  and  PcP  were  obtained  from  pro¬ 
cessing  events  with  known  magnitude  and  location.  This  was  done  for  both 
ARCESS  and  NORESS,  but  as  no  observations  were  available  for  the  FINESA 
array,  the  calibration  factors  were  temporarily  set  to  the  same  values  as  those 
of  the  ARCESS  array.  Travel  times  and  apparent  velocities  are  taken  from 
standard  tables,  and  the  receiver-source  azimuth  is  used.  The  STA  integration 
windows  are  2  seconds  for  P,  and  5  seconds  for  PP  and  PcP.  Time  tolerances 
are  2  seconds  for  P  and  3  seconds  for  PP  and  PcP. 


The  data  processing  flow  can  be  outlined  as  follows: 


•  Compute  filtered  beams  and  STA  traces  for  each  array.  The  STA  traces 
are  stored  on  cyclic  files,  containing  the  last  5  days  of  data. 

•  Compute  upper  magnitude  limits  from  the  network  of  arrays,  as  well  as 
from  each  array  separately.  The  program  reads  the  STA  traces  from  the 
cyclic  files,  and  stores  the  resulting  magnitude  traces  on  new  cyclic  files, 
also  these  containing  the  last  5  days  of  data.  The  total  size  of  the  files 
is  altogether  about  96  Mb. 

•  Analyze  the  magnitude  traces  by  an  interactive  process  within  the  IMS. 
Magnitude  data  are  read  from  the  cyclic  files  and  displayed  as  continuous 
time  series.  The  interactive  processing  is  also  attached  to  the  IMS  event 
data  base,  such  that  information  on  interesting  events  can  be  retrieved, 
e.g.  by  showing  the  location,  origin  time  and  estimated  magnitude. 

Examples  from  the  demonstration 

Fig.  7.6.2  shows  the  location  of  the  10  mining  regions  used  as  target  areas 
for  the  threshold  monitoring.  In  the  following  example,  we  will  concentrate 
on  the  mine  HC17.  In  addition  to  the  upper  magnitude  limits  derived  from 
the  network  of  arrays,  upper  magnitude  limits  were  computed  from  each  array 
separately.  Fig.  7.6.3  shows  these  traces  for  Friday  02/03/90.  The  first  6 
hours  of  the  day  are  charcterized  by  low  seismic  noise  levels,  and  any  events  at 
HC17  must  have  magnitudes  well  below  1.5.  The  low  seismic  activity  is  also 
illustrated  by  the  V’s  on  top  of  each  curve,  indicating  origin  times  of  events 
located  by  the  IMS.  An  increase  in  the  seismic  amplitude  level  caused  by  an 
event  shows  up  like  spikes  on  the  graphs.  As  expected,  these  occurrences  are 
much  more  frequent  during  working  hours  (07-15  GMT).  At  NORESS,  there  is 
also  a  general  increase  in  the  background  noise  level  during  these  hours.  This 
is  probabaly  caused  by  nearby  industrial  activity  (Fyen  et  al,  1990;  Kvasrna, 
1990). 

From  Fig.  7.6.2  we  see  that  while  the  upper  magnitude  limits  computed 
from  each  array  separately  indicate  several  time  intervals  where  events  with 
magnitudes  greater  that  2.0  may  have  occured  at  HC17,  the  network  curve 
efficiently  exclude  all  but  two  of  these  cases.  The  program  allows  us  to  expand 
the  plots,  so  that  interesting  time  intervals  can  be  investigated.  This  is  done 
in  Figure  4,  expanding  Fig.  7.6.3  for  a  time  interval  around  12  GMT.  The  two 
interesting  instances  (called  Event  1  and  Event  2)  are  identified  on  the  plot. 

On  Fig.  7.6.4  we  see  that  the  IMS  has  found  events  with  origin  times 
close  to  the  peaks.  By  clicking  on  the  V’s  with  the  mouse,  the  corresponding 
event  location  is  sent  to  the  map  for  display.  Fig.  7.6.5  shows  the  map  with 
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the  locations  of  Event  1  and  Event  2.  Event  1  is  in  fact  occnring  at  HC17, 
whereas  Event  2  is  located  approximately  400  km  further  south. 

These  two  events  illustrate  well  some  main  features  of  the  threshold  mon¬ 
itoring  technique,  and  we  will  comment  upon  them  in  some  detail: 

Event  1,  located  in  the  target  region 

In  the  case  of  Event  1,  located  at  HC17,  we  note  that  this  event  was  suf¬ 
ficiently  large  to  have  phase  detections  (both  P  and  Lg)  at  all  three  arrays. 
The  threshold  monitor  algorithm  project  the  observations  back  to  origin  times 
at  the  target  region,  so  with  an  event  occuring  at  the  target  region,  the  cor¬ 
responding  peaks  should  line  up  for  all  four  curves  in  Fig.  7.6.4.  In  addition, 
we  expect  the  upper  magnitude  limits  computed  from  each  array  separately, 
to  approach  the  same  value,  provided  that  each  array  has  peaks  exceeding  the 
background  noise  level.  For  Event  1,  both  of  these  features  are  clearly  seen.  In 
the  IMS  bulletin,  the  magnitude  of  the  event  is  estimated  to  1.6,  whereas  the 
Helsinki  magnitude  value  is  2.2.  This  difference  is  attributed  to  different  cali¬ 
bration  factors,  and  it  appears  that  the  IMS  in  general  gives  lower  magnitude 
values  than  Helsinki. 

The  network  trace  of  Fig.  7.6.4  shows  an  “upper  limit”  value  of  2.2  for 
Event  1.  It  is  important  to  be  aware  that  this  value  will  be  slightly  underes¬ 
timated  (as  a  90  per  cent  upper  limit)  in  cases  when  several  stations  actually 
detect  the  event.  The  reason  for  this  is  the  assumption  that  all  observed  phase 
magnitudes  are  less  than  or  equal  to  the  observed  “noise”  value.  In  case  of 
phase  detections  acually  attributed  to  the  event,  the  term  “less  than  or  equal 
to”  should  be  replaced  by  “equal  to”  for  those  phases,  thus  arriving  at  a  stan¬ 
dard  maximum  likelihood  magnitude  estimation  formulation  (Ringdal,  1976). 
In  practice,  this  makes  little  difference  for  a  network  of  the  type  we  are  con¬ 
sidering,  but  it  might  become  more  significant  for  larger  networks.  We  are 
currently  looking  into  ways  to  correct  for  this  bias. 

Our  conclusion  from  analyzing  Event  1  is  that  an  event  actually  occur¬ 
ring  in  the  target  region  will  have  several  readily  identifiable  features  that 
can  be  used  in  visually  confirming  the  event.  The  threshold  monitor,  when 
operated  in  conjunction  with  a  network  detector,  will  serve  to  point  out  such 
occurrences. 

Event  2,  located  outside  the  target  region 

As  seen  from  Fig.  7.6.4,  the  four  peaks  for  Event  2  do  not  line  up  very  well, 
indicating  that  the  epicenter  is  not  at  the  monitoring  region.  Note  also  that 
the  upper  magnitude  limit  inferred  from  the  network  of  arrays  is  significantly 
smaller  than  the  limits  inferred  from  each  array  separately.  From  the  map 
in  Fig.  7.6.5,  we  see  that  the  event  is  located  well  away  from  HC17.  The 
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IMS  bulletin  report  a  magnitude  of  2.6,  while  the  Helsinki  bulletin  gives  a 
magnitude  value  of  2.9.  The  upper  magnitude  limit  for  a  hypothetical  event 
at  HC17  at  the  time  considered  is  found  to  be  between  1.9  and  2.0. 

In  this  case,  the  upper  limit  represents  a  realistic  90  per  cent  confidence 
threshold  for  such  a  hypothetical  event.  Thus  we  note  that  the  method  serves 
to  ensure  that  (at  the  given  confidence  level),  a  hypothetical  event  occurring 
at  HC17  would  be  almost  a  magnitude  unit  smaller  than  the  nearby  interfering 
event. 

An  example  of  teleseismic  monitoring 

In  an  attempt  to  show  the  applicability  of  the  threshold  monitoring  method 
to  teleseismic  distances,  we  have  included  the  Shagan  River  test  site  as  a 
target  area  (Fig.  7.6.6).  The  distances  to  the  three  arrays  are  in  this  case 
between  3400  and  4300  km.  We  know  that  the  NORESS  array  has  favorable 
signal  focusing  effects  for  P-phases  from  Shagan  River,  and  that  this  will  make 
NORESS  the  most  valuable  station  for  constraining  the  magnitudes.  This 
is  clearly  illustrated  in  Fig.  7.6.7.  The  network  upper  magnitude  limits  are 
consistently  below  3.5  for  02/02/90,  and  are  only  above  3.0  when  interfering 
seismic  signals  are  observed  at  NORESS. 

We  note  that  we  do  not  have  reliable  magnitude  calibration  functions  for 
FINESA  at  the  present  time,  and  the  plots  in  Fig.  7.6.7  should  therefore  be 
interpreted  with  some  caution. 

Conclusions  and  future  perspectives 

The  implementation  of  the  threshold  monitoring  method  in  the  IMS  system 
has  shown  that  the  method  can  be  used  in  real-time  operation.  The  displays 
provided  by  the  threshold  monitor  appear  to  be  very  valuable  in  pointing 
out  time  intervals  of  particular  interest,  thus  aiding  the  analyst  in  his  work. 
The  interesting  intervals  can  be  examined  by  different  processing  techniques 
to  locate  and  identify  the  events.  Our  examples  have  demonstrated  that  the 
method  can  be  applied  both  at  regional  and  teleseismic  distances.  We  note 
however,  that  some  additional  research  needs  to  be  done  to  assess  the  potential 
bias  in  the  upper  magnitude  limits  when  detected  phases  occur  from  events  in 
the  target  region. 

A  natural  extension  of  the  implementation  would  be  to  include  more  ar¬ 
rays  or  single  stations  in  the  processing.  This  can  be  done  in  a  straightforward 
manner,  as  the  computing  algorithm  is  fully  parametrized.  Larger  geograph¬ 
ical  areas  can  be  monitored  if  standard  amplitude-distance  relationships  are 
used  to  derive  the  magnitudes.  This  will  require  careful  positioning  of  the 
target  points  and  some  research  on  the  tuning  of  the  processing  parameters. 
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Fig.  T.e.l.  Constructed  example  of  a  possible  event  interference,  where  a 
teleaeismic  P-phase  is  interfering  with  the  Pn-phase  of  a  regional  event. 


Fig.  7.6.2.  The  location  of  the  mining  regions  subjected  to  threshold  monitor¬ 
ing  is  shown  by  filled  squares.  The  three  regional  arrays,  NORESS,  ARCESS 
and  FINESA,  are  shown  by  stars.  The  filled  sectors  from  each  array  to  the 
mine  HC17,  serve  to  focus  our  interest  to  that  particular  mine  location. 
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Fig.  7.6.3.  The  top  three  panels  show  the  upper  magnitude  limits  for  mine 
HC17  for  Friday  02/02/90,  computed  from  the  three  regional  arrays  separately. 
The  unit  on  the  vertical  axes  is  magnitude  The  lower  panel  shows  the  upper 
magnitude  limit  inferred  from  the  network  of  arrays.  The  V’s  indicate  origin 
times  of  events  located  by  the  IMS. 
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Fig.  7.6.4.  This  is  a  blowup  of  Fig.  7.6.3,  with  start  time  at  11.37.31.  The 
length  of  the  time  interval  is  about  126  minutes.  The  two  marked  events  are 
of  special  interest,  since  high  upper  magnitude  limits  are  observed. 
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Fig.  7.8.5.  Locations  of  the  two  events  referred  to  in  Fig.  7.6.4. 
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Fig.  7.6.6.  This  figure  shows  the  location  of  the  Shagan  River  test  site, 
and  also  indicates  the  raypaths  (approximated  as  straight  lines)  to  NORESS, 
ARCESS  and  FINESA. 
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Fig.  7.8.7.  Upper  magnitude  limits  for  the  Shagan  River  test  site  for  Friday 
02/02/90.  The  panels  are  similar  to  those  in  Fig.  7.6.3,  and  the  V’s  indicate 
origin  times  of  events  located  by  the  IMS. 
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